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CHAPTER  I 


THE  NEGLECT  OF  QUANTITATIVE  SYNTHESIS  IN 
FEEDBACK  CONTROL  THEORY 


I  .  1  Introduction 

The  flight  control  problem  is  one  of  regulation  and  control  despite 
parameter  uncertainty  and  disturbances.  The  pilot  has  a  number  of  control 
variables  (elevator,  aileron  angles,  etc.)  available.  The  mathematical 
relations  between  these  control  inputs  and  the  output  variables  he  wishes  to 
control,  are  highly  nonlinear.  The  parameters  are  functions  of  Mach  (M)  and 
dynamic  pressure  (N)  and  are  not  precisely  known  even  if  M  and  N  are 
accurately  measured.  In  addition,  his  objectives  are  not  always  the  same. 

For  example,  in  the  longitudinal  mode  he  is  primarily  interested  In  accelera¬ 
tion  when  violently  maneuvering,  and  in  pitch  angle  when  aiming  at  a  target. 
Wind  gusts  are  external  disturbances  whose  effects  must  be  controlled. 

There  is  no  complete  synthesis  theory  as  yet  for  this  problem.  Classical 
linear  time-invariant  feedback  theory  (Nyquist  and  Bode  plots,  root  locus, 
etc.)  has  been  a  useful  tool  based,  of  course,  on  the  highly  approximate 
linear  time- invari ant  model.  Describing  function  theory  has  been  helpful  for 
the  stability  problem.  But  these  tools  are  obviously  far  from  satisfactory. 
They  must  be  accompanied  by  extensive  simulation  and  cut  and  try  modifications. 
They  have  worked  because  of  the  ingenuity  of  practical  designers  and  the 
inherent  power  of  feedback,  but  a  great  deal  of  cut  and  try  design  is  essential 
Certainly,  in  this  field,  theory  has  lagged  far  behind  practice. 
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1.2  A  Quantitative  Problem 

A  primary  reason  for  this  sad  state  is  due  to  the  almost  total  neglect 
of  a  quantitative  feedback  synthesis  theory.  After  all,  feedback  is  used  to 
achieve  desired  outputs  despite  uncertainty.  Should  not  the  extent  of  the 
uncertainty  enter  into  the  design?  For  example,  suppose  the  plant  (constrained 
part  of  the  system)  has  the  transfer  function 

kl 

P,(s)  -  -j - - -  ,  (1.1) 

s  +  A,s  +  B, 

with  uncertainties  in  k,  A,  B  given  by 

k,  e  [1,10]  ,  A]  e  [-2,4]  ,  B,  e  [2,8]  .  (1.2) 

Suppose  that  in  a  completely  different  problem,  the  plant  P2U)  has  exactly 
the  same  form  but  with  k2  £  [8,10]  ,  A2  £  [2,3]  ,  B2  e  [4,6]  .  Suppose  also 
that  the  output  performance  tolerances  are  the  same  for  both,  e.g.,  the 
system  step  response  must  lie  within  the  solid  14ne  bounds  shown  in  Fig.  1-1, 
for  all  plant  parameter  values  in  the  above  intervals.  Common  sense  suggests 
that  the  two  designs  should  be  significantly  different.  There  is  much  more 
uncertainty  in  the  first  plant.  The  "amount  of  feedback"  it  needs  should 
certainly  be  much  more  than  in  the  second  problem.  But  classical  design 
theory  ignores  this  problem.  One  presumably  emerges  with  the  same  design  for 
both  problems.  It  is  as  if  the  mere  use  of  feedback  suffices  to  scare  both 
plants  into  the  desired  behavior. 

Consider  a  variation  of  the  above  problem.  Suppose  the  uncertainties 
of  the  two  plants  are  the  same,  but  the  output  tolerances  are  different t 
e.g.,  for  one  (a)  they  are  the  dashed-line  bounds  in  Fig.  1,  and  the  solid 
ones  for  the  other  (b).  Common  sense  suggests  that  design  (a)  should  need 
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less  feedback  and  therefore  be  more  economical  In  some  sense  than  design  (b) . 
However,  again  classical  feedback  theory  cannot  cope  directly  with  this 
problem  except  by  cut  and  try,  because  it  has  no  quantitative  design  techniques. 
Modern  state-variable  design  theory  also  ignores  this  quantitative  problem. 

It  too  has  concentrated  on  design  to  achieve  a  desired  output  for  a  fixed, 
known  plant  transfer  function.  It  offers  no  design  technique  for  the  above 
problem  of  Fig.  1,  and  this  includes  the  recent  work  on  robust  design.  It 
too  must  resort  to  cut  and  try. 

This  is  indeed  a  very  sad  situation.  Feedback  is  used  primarily  to 
achieve  desired  output  tolerances  despite  uncertainty.  But  the  extent  of 
the  uncertainty  and  of  the  tolerances  do  not  appear  at  all  as  design  parameters 
in  the  vast  majority  of  the  synthesis  theories.  Thousands  of  papers  and  hun¬ 
dreds  of  books  have  been  written  on  the  subject,  but  one  hardly  finds  anywhere 
a  quantitative  problem  statement.  By  this  is  meant  a  quantitative  statement 
of  the  uncertainty  as  in  Eq.  (2),  and  a  quantitative  output  tolerance  statement 
as  in  Fig.  I.  This  makes  it  practically  impossible  to  compare  design  tech¬ 
niques,  especially  the  so-called  adaptive  techniques  which  claim  to  be  superior 
to  "conventional"  techniques. 

However,  as  described  in  Chapter  2,  quantitative  design  techniques  have 
been  developed  in  the  last  eight  years  which  can  cope  directly  with  the 
quantitative  problem  given  by  Fig.  1  and  Eqs.  (1,2).  Hence,  if  one  comes 
forth  with  claims  for  "superior"  design  techniques,  he  should  be  challenged 
to  prove  his  claims  quantitatively.  He  can  do  this  only  by  applying  his 
technique  to  a  quantitative  problem,  and  comparing  i t  to  a  design  based  on 
the  methods  of  Chapter  2.  In  fact,  the  references  noted  there  contain  many 
design  examples,  which  can  be  used  for  this  purpose.  It  is  scandalous  that 
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the  proponents  of  "adaptive"  design  and  of  "modern"  design  theories  have 
failed  to  provide  such  quantitative  proofs  of  their  claims. 

The  quantitative  techniques  described  in  Chapter  2  are  devoted  to  single 
input-output  linear  time  invariant  (It!)  uncertain  plants.  This  constitutes  a 
significant  improvement  over  conventional  nonquant i tat ive  techniques,  classical 
or  modern.  But  the  flight  control  problem  is  a  multiple  input-output  nonl inear 
problem  with  large  uncertainty,  so  the  methods  of  Chapter  2  are  also  inadequate. 
However,  there  have  recently  appeared  two  breakthroughs.  The  first  extends 
quantitative  synthesis  rigorously  (no  approximations)  to  nonlinear  uncertain 
single  input-output  systems.  The  method  is  outlined  in  Chapter  3,  and  is  applied 
in  Chapter  4  to  a  single-axis  flight  control  problem.  The  second  breakthrough 
extends  quantitative  synthesis  to  Iti  multiple  input-output  systems  with  large 
uncertainty.  It  may  be  combined  with  the  first  breakthrough,  to  permit  quanti¬ 
tative  synthesis  of  highly  uncertain  nonlinear  multiple  input-output  feedback 
systems  (Chapter  5  and  Appendix  1).  The  result  is  a  powerful  synthesis  tool, 
ideally  suited  to  the  multi -axes  nonlinear  flight  control  problem. 


Figure  1-1.  Two  sets  of  bounds  on  step  response 


CHAPTER  2 

QUANTITATIVE  FEEDBACK  SYNTHESIS  THEORY  FOR 
LINEAR  TIME  INVARIANT  SINGLE  INPUT-OUTPUT  SYSTEMS 

2.1  Origin  and  Tools 

Modern  quantitative  synthesis  theory  Is  based  to  a  significant  extent 
on  the  work  of  Bode  [I]  in  single-loop  amplifier  systems.  Bode  emphasized 
that  the  reason  for  using  feedback  was  to  control  the  plant  sensitivity 
function.  There  is  only  one  in  the  single-loop  system,  and  it  controls  the 
effect  of  parameter  uncertainty  and  external  disturbances.  Bode  showed  the 
central  role  of  the  loop  transmission  function  L(s)  and  intensively 
studied  the  properties  of  L(s)  on  the  frequency  axis,  s  =  jw  .  Unfortu¬ 
nately,  classical  feedback  control  theory  and  modern  control  theory  even  more 
so,  have  neglected  the  work  of  Bode.  His  contributions  are  still  not 
appreciated  by  most  control  specialists,  especially  that  the  essential  price 
of  feedback  is  in  the  bandwidth  of  the  loop.  The  loop  bandwidth  may  have 
to  be  hundreds  of  times  greater  than  the  control  bandwidth,  depending  on  the 
extent  of  uncertainty  and  the  tolerances  on  the  output  (recall  Chapter  1). 

In  one  extreme  (no  uncertainty)  the  loop  bandwidth  can  be  zero,  i.e.  feedback 
is  not  needed  at  all.  In  fact,  this  is  a  good  test  of  a  feedback  synthesis 
technique:  What  is  the  bandwidth  of  the  loop  transmission  when  there  is  no 
uncertainty?  If  it  is  not  zero,  then  the  design  is  not  "tuned"  to  the 
specifications  and  is  therefore  wasteful  in  loop  bandwidth.  It  is  noteworthy 
how  the  relation  of  loop  bandwidth  to  sensitivity  is  neglected  in  classical 
control  and  even  more  so  in  modern  state-variable  design  theory.  Bode  showed 
that  loop  bandwidth  reduction  is  a  central  challenge  in  feedback  theory,  and 
is  the  proper  motivation  for  nonl inear  feedback  compensation.  Thus,  an 
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adaptive  technique  should  justify  itself  by  achieving  the  same  quantitative 
tolerances  for  the  same  quantitative  uncertainty,  by  means  of  a  slower 
(smaller  bandwidth)  loop  transmission. 

2.2  Frequency  Domain  vs.  State-Domain 

It  follows  from  the  above  that  frequency  response  is  a  natural  tool  in 
linear  time- invar! ant  feedback  synthesis.  There  are  several  reasons. 

Numbers  1,  2,  6  apply  to  system  theory  as  a  whole. 

(1)  The  transfer  function  of  a  cascade  of  two  blocks  Pj  (s)  ,  P2 (s ) 
is  the  product  PjP2(s)  ,  with  their  parameters  thus  remaining  separate. 

But  the  differential  equation  for  their  cascade  connection  is  a  complicated 
mixture  of  the  individual  ones,  with  their  respective  parameters  all  mixed 
up.  If  one  wants  the  product  P ] ?2 (s )  to  have  the  properties  of  some 
desired  function  G(s)  ,  and  if  Pj  is  fixed  P2  free,  then  obviously 
P2(s)  =  G/Pj (s)  is  easy  to  determine.  This  is  not  so  if  the  cascade 
combination  is  expressed  in  differential  equation  form.  Similar  properties 
apply  to  other  kinds  of  interconnections. 

(2)  Desired  system  properties  are  relatively  easily  stated  in  the 
frequency  domain,  which  is  also  insensitive  to  the  order  of  the  system. 

They  are  much  less  easily  formulated  as  coefficients  of  a  differential 
equation,  whose  number  and  properties  are  furthermore  very  sensitive  to 
system  order. 

(3)  In  feedback  theory,  loop  bandwidth  is  the  price  paid  for  sensitivity 
reduction,  which  is  a  frequency  domain  parameter.  The  equivalent  in  the  time 
domain  is  the  speed  of  the  impulse  response  of  the  loop  transmission  and  this 
is  very  opaque  in  a  state-variable  formulation.  Integral  equations  would  be 

a  much  more  logical  time-domain  formulation  for  feedback  purposes,  but  for 
the  next  reason. 
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(1*)  In  the  problem  posed  in  Fig.  I - 1 ,  one- wants  to  control  the  response 

at  every  t  .  The  response  at  t  *  t|  is  a  function  of  all  t  <  tj  ,  the 

convolution  of  (uncertain)  functions  for  all  t  <  t |  The  problem  would 
not  be  so  bad  if  one  were  interested  in  only  a  few  tj  values,  but  as  noted 
it  is  for  all  tj  values.  In  the  frequency  domain  the  analogous  problem 
is  to  control  the  response  to  be  within  certain  bounds  (see  Section  2.3, 

Figs.  2-la,b)  for  all  u  ,  so  frequency  u>  replaces  time  t  in  the  above. 
However,  the  behavior  of  the  function  at  W|  is  very  loosely  linked  with 

its  behavior  at  all  other  o>  values.  The  constraint  of  analyticity  of  the 

s-domain  function  is  much,  much  easier  to  bear  than  the  convolution  constraint 
in  the  t-domain. 

(5)  The  nonminimum-phase  property,  as  noted  by  Bode,  is  crucial  in 
feedback  design.  It  is  explicit  and  obvious  in  the  ui-domain,  but  opaque 
and  hidden  in  the  t-domain,  whether  in  state-variables  or  in  integral 
equation  form. 

(6)  System  theory  differs  basically  from  physics,  chemistry  and  the 
specialized  engineering  disciplines  in  that  it  deals  with  the  interconnection 
of  devices  and  the  overall  system  properties  as  functions  of  the  input-output 
properties  of  the  device*.  Its  great  pride  and  boast  is  that  it  can  do 
precisely  this  (i.e.,  control  the  system's  input-output  relations),  with  no 
need  for  study  or  even  understanding  of  the  equations  and  natural  laws 
governing  the  mechanisms  of  these  individual  devices.  The  device  can  have 

a  thousand  internal  states  but  if  its  connection  to  the  overall  system  is 
via  only  a  single  input  and  a  single  output,  then  the  system  theorist  needs 
only  the  single  equation  relating  these  two  variables.  What  is  the  point  of 
cluttering  up  the  formulation  with  the  thousand  internal  variables?  The 
essence  of  good  engineering  is  simplicity  and  economy  of  representation,  In 
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order  that  one  may  see  the  forest  from  the  trees.  The  problems  are  usually 
difficult  enough  without  this  unnecessary  clutter. 

(7)  It  has  been  argued  that  the  state-variable  description  is  more 
accurate  and  fundamental.  Obviously,  it  is  more  detailed,  but  why  does  this 
make  it  more  fundamental?  The  latter  is  a  subjective  concept  related  to 
one's  objectives.  As  George  Zames  has  noted,  that  in  pursuing  this  argument 
a  quantum  mechanical  description  is  far  more  accurate  and  "fundamental"  than 
state-variable. 

(8)  Much  has  been  made  of  the  concepts  of  "controllability  and 
observability"  as  a  justification  for  the  state-variable  description  and 
methodology.  But  it  has  been  shown  [2]  that  a  system  can  be  controllable  and 
observable,  and  yet  is  totally  inadequate  for  proper  control.  Thus,  these 
concepts  can  be  highly  misleading.  Also,  it  has  been  shown  they  are  not  at 
all  needed  even  for  what  they  purport  to  do,  if  one  formulated  the  problem 
with  the  concept  of  uncertainty. 

(9)  Much  has  been  made  of  the  fact  that  the  Laplace  transform  is 
unavailable  in  linear  time-varying  and  in  nonlinear  systems,  whereas  the 
state-variable  formulation  is  universal.  To  this  the  answer  is,  "So  what?" 
The  objective  is  synthesis,  not  mere  formulation.  Quantitative  techniques 
with  deep  insight,  for  the  synthesis  of  uncertain  systems  to  satisfy  quanti¬ 
tative  performance  specifications,  are  available  in  the  transform  domain, 
for  linear  time-invariant  systems.  They  are  not  available  by  state-variable 
methods.  Of  what  value  is  a  universal  formulation  if  it  cannot  cope  with 
the  quantitative  uncertainty  problem  even  in  its  simplest  category — the 
linear,  time- invariant  one. 


Ironically,  transform  concepts  have  led  to  precise  quantitative  design 
techniques  for  both  time-varying  and  nonlinear  systems  [3~5l •  The  application 
of  such  a  technique  to  the  flight  control  problem  is,  in  fact,  the  main 
theme  of  this  report. 

It  should  be  emphasized  that  there  are  areas  in  system  theory  where  the 
state-variable  formulation  is  the  natural  one,  e.g.  in  optimal  control  if 
the  realistic  cost  function  is  indeed  a  function  of  all  the  states.  It  is 
not  our  purpose  to  derogate  state-variables,  but  rather  to  present  honestly 
the  role  of  frequency  response  in  quantitative  feedback  theory.  In  such  a 
presentation  one  naturally  contrasts  it  with  stale  variable  theory.  At  least 
our  presentation  gives  detailed  arguments,  whereas  frequency- response  is 
generally  cast  aside  by  state-variable  enthusiasts,  with  hardly  any  detailed, 
reasoned  discussion.  Also,  our  argument  is  not  against  state-variable 
Formulation  per  se,  but  to  the  synthesis  technique  that  is  used  thereafter. 

2.3  Translation  of  Time  Domain  into  Frequency  Domain  Specifications 

In  the  synthesis  techniques  listed  next,  there  is  always  a  constrained 
Plant  which  is  described  by  a  system  of  differential  equations  whose  para¬ 
meters  are  uncertain,  giving  a  set  of  plant  functions  P  =  { P }  .  Thus  in 
Eqs.  (1.1,  2),  each  parameter  combination  gives  a  different  P  e  P  .  The 
same  design  philosophy  can  be  used  for  sampled-data  systems  [6], 

The  objective  is  to  achieve  certain  apriori  specified  performance 
objectives  V  P  e  P  .  If  the  overall  system  is  to  be  linear  time- invariant 
(lti)  even  if  the  plant  itself  is  nonlinear,  it  can  be  characterized  by  its 
response  to  any  input,  and  the  step  response  is  very  popular  because  it 
combines  within  it  both  the  fastest  kind  of  input  (an  abrupt  change)  and  the 
slowest  (no  change).  Time  domain  specifications  are  reasonable  in  many  cases, 
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as  in  Figure  2-la,  where  the  step  response  is  to  be  inside  the  bounds 
bj,b2  VP  e  P  ,  with  additional  bounds  of  similar  nature  on  the  first  and 
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Figure  2-la.  Time  domain  step  response  specification. 


perhaps  higher  derivatives.  Our  design  technique  is  in  the  frequency-domain, 
so  we  must  translate  such  t-domain  bounds  into  "equivalent"  w-domain  bounds 
on  the  system  frequency  response  T(ju))  .  If  the  system  is  minimum-phase 
[6],  j T ( joj) |  suffices  and  we  restrict  ourselves  here  to  such  systems. 

This  translation  is,  as  of  this  date,  an  engineering  art  rather  than  a 
science.  Advice  on  how  to  translate  is  scattered  in  the  literature  [ 7_ 9 ] - 
Very  good  results  have  been  obtained  with  only  moderate  effort.  We  shall 
assume  in  this  work  that  the  translation  has  already  been  done.  It  Is  worth 
noting  that  it  has  been  shown  [4]  that  for  minimum-phase  systems,  time-domain 
specifications  on  the  step  response  and  on  its  derivatives  of  the  following 
nature 

b2(0(t)  <  c(i)(t)  <  bj^U)  ,  i  -  0,1 . n,te[0,-)  (2.1) 


-10- 


•JO.  9 
<*  -  • 


-20 


^Vjx  B, 

8,'db\ 

\ 

B2  \ 


Figure  2-lb.  Frequency  domain  specification. 

can  always  be  satisfied  by  means  of  urdomain  bounds  of  the  following  nature 

B2(w)  £  } T ( joj)  |  _<  Bj  (oj)  .  (2. 

Hence,  it  is  guaranteed  that  there  exist  w-domain  bounds  of  the  form  of 
Fig.  2-lb,  which  satisfy  time-domain  bounds  of  the  form  Eq.  (2.1). 

An  example  of  time-domain  bounds  and  their  equivalent  u)-domain  bounds 
is  shown  in  Figs.  2-!c,d.  Fig.  2-lc  also  shows  the  simulation  results 
obtained  from  the  design,  which  was  incidentally  a  multiple-loop  one  (see 
Section  2.7),  for  the  structure  of  Fig  2-1 3b .  Such  good  correlation  between 
t  and  u)  domain  bounds  are  not  atypical. 
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2.4.  The  single-loop  two-degree-of- freedom  system  [7].  (Fig.  2-2) 

Here  the  plant  output  c(t)  is  assumed  measurable  and  available  for 
feedback  and  so  is  the  command  input  r(t)  .  The  processing  of  these  two 
signals  provides  two  independent  compensation  functions  to  the  designer. 

An  infinitude  of  canonical  two-degree-of- freedom  structure  may  be  used  [6]. 
The  design  procedure  developed  in  [7]  used  Figure  2-2a,  but  suppose  the 


0 
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Figure  2-2.  Structures  of  2-D.O.F.  system. 


sensor  transfer  function  is  H(s)  ,  then  one  can  use  Figure  2- 2b ;  letting 
GjH  (of  Figure  2- 2b ) )  =  G  (of  Figure  2-2a)),  in  order  to  have  the  same 
loop  transmission  function  l(s)  *  GP  *  G j PH  ,  and  F j G ^  =  FG  in  order  to 
have  the  same  system  transfer  function 


(s)  _ 

RlsT  ’ 


T(s) 


FGP 
I  +  GP 


CCD 

Vr 

1  +  g  ,  PH 


(2.3) 
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2.4.1  A  2-degree-of- freedom  structure  with  2-loop  implementation 


Suppose  large  loop  feedback  bandwidth  is  needed  and  it  is  found  that  an 
independent  sensor  measuring  c(t)  (e.g.,  a  tachometer  in  a  position  servo) 
gives  less  noise  than  the  differentiation  of  a  position  sensor,  so  both 
sensors  are  used,  as  in  Figure  2-3,  with  the  two  sensor  transfer  functions 
Hj  ,  H2  ,  and  say  the  structure  in  Figure  2-3  is  used.  This  is  a  two-loop 


l. 


Figure  2~3-  2-loop,  2-D.O.F.  structure. 

structure  physically,  but  in  terms  of  fundamental  feedback  design  it  is  a 
two-degree-of-f reedom  system,  so  the  quantitative  design  theory  of  Figure  2-2a 
is  used,  giving  G  and  F  .  It  is  required  that  the  loop  transmission  around 
P  ,  be  the  same  in  both  cases,  i.e., 


and 


so 


L  =  PG  (Fig.  2-2a)  =  P  [G  (H,  +  GbH2) ] 


T  =  F  I  +  GP  (Fi9,  Z'2a}  =  Fb  I  +PlGa(H,a+  GbH2)J 


(2.4a) 


(2.4b) 


G  =  Ga(H)+GbH2)  ,  FG  =  F^  . 


(2.4c) 


Hj  ,  H2  are  known,  so  one  must  decide  how  to  spilt  G  ■  Ga (H |  + G^Hj)  between 


Ga  and  Gb  .  This  is  done  by  considering  the  effect  of  sensor  noise  Nj  ,  N2 
at  the  plant  input. 


given  that 


-XN(ja)) 


Ga(HlVH2SbN2) 

1  +P[Ga(H,  +GbH2)] 


(2.5) 


G  =  £  (jw)  -  Ga(H,  +H2Gb) 


(2.6) 


is  fixed  by  the  quantitative  design  technique  of  [7l • 


The  objective  is  to  minimize 


subject  to  the  above  con¬ 


straint.  This  is  a  straightforward  optimization  problem  which  can  be  solved 
outside  the  realm  of  quantitative  feedback  synthesis.  The  latter  only  provides 
the  design  with  the  feedback  loop  transmission  (L)  needed  around  the  plant, 
and  the  prefilter  (F)  needed  to  process  the  command  input  r(t).  The  state- 
of-the-art  in  sensors  and  in  filter  synthesis  determines  how  L  and  F  are 
to  be  realized.  In  fact,  in  the  above  context  one  might  consider  use  of  an 
accelerometer  in  a  3-loop  feedback  structure.  But  from  our  point  of  view 
the  structure  remains  that  of  a  two-degree-of- freedom  system  and  we  shall 
continue  to  associate  the  latter  with  a  single-loop  system.  The  single-loop 
design  technique  is  the  basic  building  block  for  ail  the  other  more  complex 
structures,  so  it  is  next  reviewed. 


2.1*. 2  Review  of  Two- Deg ree-of- Freedom  Quantitative  Design  Theory 

GP 

Figure  2-2a  is  used  with  T  »  F  y+'Qp’  •  It  is  assumed  that  the  compen¬ 
sation  networks,  whose  power  levels  can  be  very  low  (as  the  plant  contains 
the  power  elements),  can  be  constructed  with  negligible  uncertainty  in  their 
transfer  functions.  Hence,  due  to  the  uncertainty  in  P  , 


fm  1.1 


and 


A  Jin  T  =  A  Jin  =  A  Hn  •—-  ,  L  -  GP 


(2.7) 


A  Jin  |  T  ( jto)  J  »  A  Jin 


L  (jui) 

r+  ujwj 


(2.8) 


Given  that  the  specifications  require  that  A£n | T ( jto)  |  _<  6jdb  at  u> j  in 
Figure  2-lb,  what  are  the  resulting  constraints  on  L(jujj)  7  It  is  convenient 
to  pick  a  "nominal"  plant  Pq ( s )  ,  and  derive  the  bounds  on  the  resulting 
"nominal"  loop  function  Lg  =  PgG  .  These  bounds  can  be  found  by  means  of 
a  digital  computer,  but  it  is  very  useful  for  insight  to  see  it  done  on  the 
Nichols  chart  (logarithmic  complex  plane  with  abscissa  in  degrees,  ordinate 
in  decibels  =  2C  log^).  The  procedure  is  illustrated  for  the  case 

F ( s )  —  y  ;  ke  [1,10],  a  e  [1,10]  ,  (2.9) 

and  say  k  =  1  ,  a  =  10  are  chosen  as  nominal,  giving  PQ  =  10/(s(s  +  10)). 

At  w  =  2  rps,  F  ( j  2 )  lies  within  the  boundaries  given  by  AECD  in  Figure  2-4. 
Since  Jin  L  =  Hr,  G  +  £n  P  ,  the  pattern  outlined  by  ABCD  may  be  translated, 
but  not  rotated,  on  the  Nichols'  chait,  the  amount  of  translation  being  giver, 
by  the  value  of  Jin  G ( j 2 )  .  For  example,  if  a  trial  design  of  L ( j 2 ) 
corresponds  to  the  template  P ( j 2 )  at  A'B'C'D1  in  Figure  2-4,  then 

Glj2)  db  -  IUJ2)ldt  -  |P(J2)ldb 

U.  10/ 

=  (-2.0)  -  (-13.0)  =  11.0  db 


Arg  G ( j 2 )  =  Arg  L(j2)  -  Arg  P(j2) 

«  (-60°)  -  (-153.4°)  =  93.4°  . 


(2.11) 
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so  the  above  trial  position  of  | Lq ( j 2) |  is  in  this  case  more 
than  satisfactory.  The  template  is  lowered  on  the  Nichols'  chart  to  A"B"C"D", 
where  the  extreme  values  of  An|L/(l+L) |  are  at  C"  (-0.7  db)  ,  A"  (-7-2  db). 
Thus,  if  Arg  Lg(j2)  *  -60°  ,  then  -4.2  db  is  the  smallest  magnitude  of 
Lq ( j 2)  which  satisfies  the  6.5  db  specification  for  A  An  |T|  .  Any  larger 
magnitude  is  satisfactory  but  represents  over-design  at  that  frequency. 

The  manipulation  of  the  u  =  2  template  is  repeated  along  a  new  vertical 
line,  and  a  corresponding  new  minimum  of  | Lq ( j 2 ) |  found.  Sufficient  points 

are  obtained  in  this  manner  to  permit  drawing  a  continuous  curve  of  the  bound 
on  Lq ( j 2)  ,  as  shown  in  Figure  2-4.  The  above  is  repeated  at  other  frequencies, 
resulting  in  a  family  of  boundaries  of  permissible  Lq(Jo))  . 

2.4.4  Nature  of  the  bounds  on  LqQcj) 

A  typical  set  of  bounds  is  shown  in  Figure  2-5-  The  bounds  tend  to 
move  down  in  the  Nichols  chart  (become  less  onerous),  obviously  because  as 
a)  increases,  greater  change  in  |T(jw)|  is  permitted,  as  in  Figure  2-lb. 
However,  they  may  cross  and  even  temporarily  move  up  higher  with  increase 
of  <jj  .  It  is  in  fact  essential  that  at  large  enough  io  ,  the  uncertainty 
in  |T(ju))|  (i-e.,  the  bounds  on  |  T  ( joo)  |  )  be  greater  than  the  uncertainty 

in  P (joo)  ,  because  the  net  sensitivity  reduction  is  always  zero  in  any 
practical  system  as  was  long  ago  [1]  shown  by  Bode, 

/»  oo  -  00 

An  |  si ( jco)  |  du  =  -  An  1 1  +  L(ju>)  |  dw  =  0  (2.12) 

J0  J0 

where  Sp  =  9pyp  is  the  sensitivity  function. 

L  a 

In  the  above  examples  as  w  -*■  ®  ,  P  -►  -y  ,  so  A  An  |  P  |  -*■  A  An(ka)  ■  40  db. 

s4 

Note  in  Figure  2-lb  that  the  permitted  A  AnjT(jeo)|  »  40  db  as  u>  >  50  . 
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Figure  2~5.  Typical  bounds  on  Lg  in  Nichols  chart. 

Such  large  tolerances  on  J T ( jo>)  j  at  large  u>  are  tolerable  because 
|T(juj)|  is  negligible  at  large  w  ,  e.g.,  if  |P(ju>)|  can  change  at  most 

by  kO  db  at  large  u)  but  J T ( jo>)  |  changes  by  52  db  ,  who  cares  if  this 

52  db  change  is  from  1 T |  ■  10  ^  to  |Tl  w  *  400  x  10  ^  .  In  return, 

one  can  concentrate  the  sensitivity  reduction  over  the  bandwidth  of  T(jw)  . 
Thus,  although  |P(ju))J  in  this  region  varies  by  say  40  db  ,  | T ( ju>)  |  may 

be  controlled  to  vary  by  only  U  db  ,  or  0.04  db  if  desired. 
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2.4.5  Universal  high-frequency  (UHF)  boundary 

As  noted,  in  the  high-frequency  range  A  2n|T(juj)|  must  realistically 
be  allowed  to  be  >>  A  £n|P(ju)j  ,  and  this  is  reflected  in  the  bounds  on 
L0{jw)  tending  to  a  very  narrow  pencil  B^  in  the  arithmetic  complex  plane 
*Sni 

P0  -*■  — “  as  s  ->  »  )  as  in  Figure  2-6a  and  as  in  2-6b  in  the  Nichols' 
s 

chart.  In  Figure  2-6b,  the  boundary  B.P  is  drawn  for  the  case 


Figure  2.6a.  Typical  bounds  on  L0(jw)  in  complex  plane. 


A  Jin  L  3  A  In  k  »  20  db  ,  A  fcn|T(jw)  |  =  A  2n|L/(l+L)  |  <  23  db  at  . 
However,  the  resulting  peak  value  of  |L/(1+L)|  is  23  db  -  IA.1  arithmetic 
^  "  ^max  ’  indicating  a  highly  under-damped  pole  pair  at  the  corresponding 


frequency  with  damping  ratio  K  *  0.03**  >  when  k  *  kmax  .  This  tremendous 
peaking  does  not  appear  in  the  system  response  to  the  command  inputs  R  , 
because  it  is  filtered  out  by  the  pre-filter  F  in  Figure  2-2a.  But  the 
system  response  to  a  disturbance  D  in  Figure  2-2a,  is  given  by 
Td  *  TJ  “  (*+L)  •  Disturbance  attenuation  generates  its  own  requirements 

on  L  ,  which  may  lead  to  more  stringent  bounds  on  L  than  those  due  to 
T(j  )  .  The  final  contours  used  in  the  design  [7]  must  be  the  most  stringent 
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Figure  2-6b.  Bounds  on  Lg(jw)  on  Nichols'  chart 
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composite  of  the  two.  However,  even  if  D  is  very  small,  it  is  usually 


certain  that  a  peak  |T^|  of  14.1  is  intolerable.  It  is  reasonable  to  add 
a  requirement  |  |  £  Y  a  constant,  for  all  u)  and  over  the  whole 

range  of  P  parameter  values.  The  resulting  constraining  contours  denoted 
by  are  shown  in  Figure  2-6b  for  the  case  A  Jin  k  ■  20  db  ,  and  for 

Y  =  2. 3  ,  3- 5  ,  5  db  (all  those  contours  are  symmetrical  with  respect  to  the 
vertical  Arg  L  ■  -180°  on  the  Nichols'  chart).  I f  y  ■  5  db  is  used, 
then  B{u>v)  indicates  the  composite  contour  shown  in  Figure  2-6b.  For 
w  >  a)  ,  |A  T(ju>)|  increases  while  y  remains  the  same,  so  that  sooner  or 
later  there  is  reached  a  frequency  3  B(oj)  =  Bh(y)  ,  v  w  il  •  This 
boundary  is  called  the  "universal  high  frequency"  (UHF)  boundary. 

2.4.6  Th~  optimum  L ( juo) 

It  has  been  shown  [10, 11]  that  a  realistic  definition  of  optimum  in 

the  lti  system  is  the  minimization  of  k,  defined  by  limL(s)  =  ks  e  , 

s  oo 

where  e  is  the  excess  of  poles  over  zeros  assigned  to  L(s)  . 

It  has  been  proven  [10,11]  that  the  optimum  L  lies  on  its  boundary 
B.  at  each  w.  and  that  such  an  optimum  exists  and  is  unique.  Most 
important  for  the  present  purpose,  is  that  in  significant  plant  ignorance 
problems  the  ideal  optimal  L  has  the  properties  shown  in  Figure  2-7,  i.e., 
over  a  significant  range  it  follows  along  UV  up  to  the  point  J  at 

which  it  abruptly  jumps  to  infinity  along  WW'W"  and  returns  on  the  vertical 
line  YZ  ,  whose  phase  is  (-90°)  e  .  Such  an  ideal  L(jio)  is,  of  course, 
impractical.  A  practical  suboptimum  L  is  shown  in  Figure  2-7. 

Some  results  of  a  numerical  design  example  are  shown  in  Figure  2-8. 

They  were  derived  for  the  following  problem. 
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Figure  2-7.  Bounds  on  Lg  and  Optimum  Lg  on  Nichols’  Chart. 

2.4.7.  Numerical  example.  (Figures  2-2a,  2-8,  Section  2.6) 

Plant:  P  -  P^  -  Pj  (P£Pp 

Kl 

Plant  uncertainty:  P.  -  -5 - =•  ,  "3  <  £10  _<  5 

1  S*  +  2&)  S  +  u  n 

n  n 

2  <  U  <  10  ,  4  <  K.  <  1250 

—  n  —  —  1  — 

I  <  A  <  3  ,  10  <  K£  <  33-3 

10  <  B  <  20  ,  100  <  <  158 

Performance  Specification:  Shown  in  Figure  2-lb,  were  originally  derived 

from  time  domain  bounds  [9]. 


P ' 

K3 


AK£ 

S+A 

BKj 

S+B 


Disturbance  response: 


y  <  3-0  db 


n 


-3oc’  -240°  >180°  _ - 120** 


Figure  2-8.  Single  loop  LsC  and  outer  loop  L^Q 
of  the  numerical  example  in  Section  2.^.7- 
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This  example  is  used  in  Section  2.6  in  a  demonstration  of  the  advantages  of 
multiple-loop  design. 

2.5  Cost  of  Feedback  and  Effect  of  Sensor  Noise 

In  significant  plant  ignorance  problems,  there  is  a  strong  tendency 
for  the  design  to  be  such  that  N  ,  in  Figure  2-2a,  is  so  highly  amplified 
as  to  saturate  the  plant  input  at  X  .  The  noise  response  function  is  (see 
Figure  2-2a) 

T  £  X  _  -C  -L/P 
N~N~  1  +  GP  =  1+L 

=  -L/P  in  h.f.  range  where  |L|  «  1  .  (2.13) 

N  represents  the  square  root  of  the  noise  power  spectrum.  The  noise  response 
of  the  numerical  design  example  of  the  last  section  is  shown  in  Figure  2-3* 
Notice  that  the  noise  component  at  x,  in  Figure  2- 2a,  is  most  important 
in  the  high-frequency  range  where  the  useful  command  and  disturbance  components, 
due  to  R  and  C  are  relatively  small,  rather  than  in  the  low  frequent^ 
range  where  the  latter  are  relatively  large.  Hence,  it  is  desirable  tc 
decrease  |l|  vs.  w  ,  as  fast  as  possible  in  the  high  f requency  range. 

Ever,  a  saving  which  is  small  in  the  logarithmic  scale  car.  be  significant 
in  rms  sensor  noise  effect.  Reduction  in  this  'cost  cf  feedback'  is  the 
primary  motivation  for  turning  tc  multiple-loop  design. 
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Figure  2-9.  Enormous  sensor  noise  amplification  in 
the  single-loop  design  for  Section  2.4.7  example. 

2.5. 1.  Reduction  in  Cost  of  Feedback  by  means  of  linear  time-varying 
compensation  and  nonlinear  compensation. 

To  reduce  the  hf  sensor  noise  effect,  one  way  is  by  linear  time- 
varying  compensation  if  the  problem  has  time-varying  features  [12].  Another 
is  by  nonlinear  compensation.  Actually  the  so-called  "adaptive"  system  is 
in  the  category  of  nonlinear  compensation.  They  may  or  may  not  be  better 
than  lti  compensation  in  reducing  the  'cost  of  feedback'.  It  is  noteworthy 
and  scandalous  that  in  the  vast  literature  on  adaptive  systems,  there  is 
hardly  ever  any  quantitative  comparison  between  the  adaptive  design  promoted 
and  a  proper  lti  design  accomplishing  the  same  design  objectives.  One  could 
excuse  this  not  being  done  in  a  general  manner,  because  there  is  hardly  any 
'adaptive'  method  permitting  quantitative  design  in  the  sense  here  defined. 
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However,  it  could  at  least  be  done  experimentally.  Occasionally  one  sees  a 
comparison,  with  an  'ordinary'  or  so-called  'classical'  design.  But  the 
comparison  is  usually  greatly  biased,  because  generally  some  very  naive  lti 
design  is  used,  and  there  is  no  statement  of  specifications — even  made  up 
after  the  fact.  There  is  not  recalled  a  single  comparison,  on  the  part  of 
the  proponents  of  adaptive  systems,  with  the  lti  quantitative  design  tech¬ 
nique  [7]  discussed  here.  Some  nonlinear  compensation  techniques  for  which 
a  quantitative  design  theory  exists  to  a  greater  or  lesser  extent  have 
appeared  in  the  literature  [13-17l  for  which  such  comparisons  are  possible. 

It  is  noteworthy  that  these  were  expressly  motivated  by  the  desire  to  reduce 
the  'cost  of  feedback1,  so  that  such  comparisons  were  a  natural  by-product. 

2.6  Mui t iple- loop  feedback. 

Another  method  of  'cost  of  feedback'  reduction,  in  the  context  of  lti 
design,  is  by  means  of  multiple-loop  feedback,  restricted  to  those  cases 
where  an  additional  plant  variable  (besides  the  plant  output)  are  available 
for  feedback  purposes.  Such  a  multiple  loop  design  technique  was  first 
developed  [ 1 3]  for  the  cascaded  structure  of  Figure  2-10a. 

The  design  example  of  Section  2.4.7  was  done  by  means  of  a  two-loop 
design  (n  =  2  in  Fig. 2-10).  The  resulting  new  outer  loop  (Ljq)  is  shown 
in  Figures  2-8,  2-11  and  is  considerably  more  economical  in  bandwidth  than 
the  single-loop  design  (L^g)  for  the  same  problem.  The  effects  of  the  outer 
sensor  noise  Nj  at  the  plant  input  are  compared  in  Fig.  2-12.  Note  the 
tremendous  improvement.  However,  there  is  now  a  second  sensor  with  noise 
source  N£  in  Fig.  2-1  Ob.  Its  effect  is  shown  in  Fig.  2-12.  It  can  be 
reduced  by  using  3  loops  (n  =  3  in  Fig.  2-10).  Compare  X/N^  for  2  and 
3  loop  designs  in  Fig.  2-12.  But  now  there  is  Nj  to  consider  and  X/N^ 
is  shown  in  Fig.  2-12. 
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2.6.1.  Design  perspective 

The  above  results  are  very  impressive.  But  what  is  even  more  impressive 
is  that  a  very  good  approximation  to  the  final  mul tipi e- loop  designs  and  the 
sensor  noise  effects,  can  be  gotten  very  quickly  by  means  of  a  technique 
called  "Design  Perspective"  [19,20].  The  application  of  Design  Perspective 
to  the  above  example  is  shown  in  Fig.  2-11.  A  detailed  design  was  required 
only  for  the  single-loop  L^q  .  The  dashed  lines  show  the  results  using 
Design  Perspective  for  Ljq- outer  loop  in  2  or  3~loop  design,  l^'first 
inner  loop  in  a  2-loop  design,  -first  inner  loop  in  a  3-loop  design, 

L^- second  inner  loop  in  a  3-loop  design.  Note  the  excellent  agreement 
between  the  approximate  Design  Perspective  results  and  the  final  detailed 
designs.  Design  Perspective  enables  the  designer  to  obtain  a  good  under¬ 
standing  and  perspective  of  the  important  design  trade-offs  quite  early  in 
the  game,  without  the  need  of  a  detailed  design. 

2.7  More  Complex  Multiple-Loop  Structures  Including  Plant  Modification 

Systematic  quantitative  design  techniques  have  also  been  developed 
[19,20]  for  the  more  complex  mul tiple- loop  structures  of  Figures  2-1 3a, b 
denoted  as  the  Triangular  and  Parallel-Cascade  structures  respectively. 

The  results  shown  in  Figures  2-lc,d  were  obtained  for  a  design  example  based 
on  Figure  2-1 3b  with  two  parallel  sections  (m  =  2)  nj  =  3  ,  n2  =  1  and 
Pj.  *  kj./s  f°r  '  =  1,2,3,  P21  =  k2/s^  ,  and  very  large  uncertainties 

kn  c  [50,5003  ,  k]2  e  [20,800]  ,  k]3  e  [1,60]  ,  k£  e  [1000^,200,000]. 

The  parameter  values  in  the  brackets  in  Fig.^-lc  have  been  normalized  and 
are  given  in  the  order  k^  .  |  »  ^13  »  ^2  ’  Design  Perspective  has  also 

been  developed  for  these  structures  [19,20]. 
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Plant  Modi fi cat ion 


It  is  noted  that  in  Figures  2-10,  I3a,b,  sach  feedback  loop  is  returned 
to  the  plant  input  X  .  No  feedback  is  allowed  to  any  internal  plant 
variables  e.g.  from  C  to  C ^  or  more  generally  from  any  C.  to  Cj  , 
j  >  i  and  j  =  2,3,...,n-  I  .  The  reason  is  that  any  such  internal  feedback 
constitutes  "plant  modification".  The  plant  has  been  assembled  by  its 
specialists  to  deliver  some  maximum  output  C  and  the  permissible  levels  of 
C2  =  C/Pj  ,  =  C/P j ?2» • • • ,C.  =  C/Pj  ...  P._j  are  thereby  determined. 

Suppose  there  is  feedback  from  C  vis  Hj  to  Cj  ,  as  shown  in  the  insert 

r 

in  Figure  2-10a.  Now  X2  =  p—  (1+PjHj)  with  signal  level  possibly  much 

greater  than  the  previous  C/P  j  ,  which  the  plant  may  perhaps  not  be  able 
to  handle. 

We  thus  assume  that  the  "feedback  specialist"  is  called  in  to  design 
the  feedback  network  around  the  plant,  after  the  latter  has  been  built  by 
the  "plant  specialist".  This  is  the  situation  very  often  in  practice.  If 

the  feedback  specialist  does  his  job  properly,  i.e.  achieves  the  system 

response  function  T(s)  within  its  tolerances  V  P  £  P  ,  then  the  signal 

levels  inside  the  plant  will  be  within  the  values  allowed  by  the  plant 

specialist,  so  long  as  the  command  input  functions  r(t)  are  in  the  set  for 
which  the  system  was  designed. 

Recently,  quantitative  design  has  been  extended  for  the  first  time 
[21-2],  specifically  to  the  single  branch  cascade  plant,  with  plant  modifica¬ 
tion  allowed.  Some  of  the  resulting  structures  are  shown  in  Figure  2-l,»- 
The  degree  of  modification,  in  a  rms  sense,  of  the  internal  plant  variables 
was  added  as  one  more  design  specification,  in  addition  to  those  listed  here. 
It  was  shown  that  the  loop  bandwidths  can  thereby  be  significantly  decreased, 


[ure  2-1 Some  feedback  structures  with  Plant  Modification. 


beyond  that  possible  in  "no  plant  modification"  designs.  This  indicates 
that  in  significant  plant  uncertainty  problems,  it  is  definitely  advantageous 
to  have  the  "feedback  specialist"  participate  with  the  "plant  specialist" 
in  the  design  of  the  plant. 


2.8  Summary 

Some  of  the  principal  features  of  single  input-output  linear  time 
invariant  ( 1 1 i )  quantitative  synthesis  theory  for  highly  uncertain  plants, 
have  been  presented  briefly  in  this  chapter.  It  should  be  clear  that  a 
respectable,  mature  synthesis  theory  exists  for  this  class,  although  much 
remains  to  be  done.  Numerous  quantitative  design  examples  have  been  executed 
together  with  computer  simulations,  which  have  corroborated  the  design 
techniques.  The  existence  of  these  design  examples  serve  as  benchmarks 
against  which  any  "adaptive"  design  with  claims  of  superiority  can  be  checked. 
These  design  techniques  can,  of  course,  be  applied  to  the  lti  flight  control 
problem.  But  what  is  of  much  greater  importance  is  the  startling  fact  that 
most  of  these  techniques  are  applicable  exactly  (no  approximations)  to  linear 
time-varying,  nonlinear  and  even  nonlinear  time-varying  highly  uncertain 
plant  problems.  This  is  treated  in  the  next  chapter. 
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CHAPTER  3 

A  BREAKTHROUGH  IN  QUANTITATIVE  FEEDBACK  SYNTHESIS 
3.1  Introduction 

Quantitative  feedback  synthesis  in  the  frequency  domain  has  clearly 
shown  its  value  in  the  single  input-output  multiple-loop  design  techniques 
for  lti  systems,  which  were  briefly  outlined  in  Chapter  2.  However,  a 
fantastic,  hitherto  considered  impossible  achievement  has  recently  been 
made  [3~5] -  It  has  been  rigorously  proven  that  these  techniques  can  be 
applied  exactly  to  a  large  class  of  nonlinear,  highly  uncertain  plants,  even 
nonlinear  time- varying  plants.  The  interesting  point  is  that  factors  like 
uncertainty,  synthesis,  quantitative  design  to  specifications  are  generally 
considered  to  make  the  design  problem  much  more  difficult.  But  it  was  pre¬ 
cisely  the  precision  and  discipline  involved  in  such  strict  requirements, 
which  led  to  the  development  of  the  rigorous  nonlinear  synthesis  technique. 
This  claim  of  precise  design  of  nonlinear,  uncertain  systems  seems  so 
unbelievable  that  it  is  worth  repeating  that  there  are  no  approximations 
involved,  no  linearizations,  no  describing  function  type  of  approximations, 
etc. 


3.2  Concept  of  the  Linear  Time  Invariant  Equivalent  Set  (LTIE) 

The  simple  single-loop  structure  of  Figure  3"la  is  first  treated,  for 
the  case  of  a  single  system  command  input  r^(t)  .  The  nonlinear  plant  w, 
with  y ( t)  =w(x(t))  ,  has  uncertain  parameters,  e.g.  w  is  given  by 

y  +  Ay^y  +  Bym  s  i  gn  y  =  Ex  (3.1) 

with  A,  B,  m,  E  uncertain  in  that  it  is  only  known  that  A  e  [-1,3]  , 
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Figure  3-la.  Single  loop  feedback  system  with  nonlinear  uncertain  plant  w  e  W. 
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Figure  3~lb.  Tolerances  on  response  to  2-unit  step  command. 


B  e  [2,8]  ,  E  e  [1,2]  ,  m  e  [.3,1.2]  .  Each  combination  of  possible  values 
of  A,...,m  give  a  different  w,  thus  generating  a  set  W=  {w}  of 
nonlinear  plant  functions.  It  is  required  that  in  response  to  command  input 
r j ( t )  ,  the  system  output  y(t)  should  be  a  member  of  a  specified  acceptable 
set  A  .  For  example,  if  r ^ ( t )  is  a  step  function  of  2  units  magnitude, 
the  tolerances  on  the  output  and  its  derivative  may  be  as  shown  in  Fig.  3“lb. 
Any  y(t)  which  satisfies  these  tolerances  is  a  member  of  the  acceptable 
set  A  *  {a  ( t ) }  . 

The  next  step  is  to  find  a  linear  t ime- invariant  plant  set  P  ,  which 
is  precisely  equivalent  to  the  nonlinear  plant  set  W  ,  for  the  purpose  of 
the  synthesis  problem.  To  do  this,  take  any  a.  e  A  and  any  plant  w-*  e  W  . 
Find  the  input  to  wJ  which  gives  the  output  a.  ,  i.e.  solve  Eq.  (3-1)  for 
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x  given  y  ,  which  in  this  case  is  a  relatively  easy  operation.  Denote  the 

resulting  x  by  xj  .  Now  find  a  linear  time-invariant  plant  P-j  ,  which 
is  equivalent  to  w*  for  this  special  case  i.e.  the  output  of  P-j  is  a.(t)  , 
when  its  input  is  x-!(t)  .  An  easy  way  to  do  this  is  to  let 


This  implies  that  both  a.  (t)  and  x-j  (t)  have  Laplace  transforms,  a  condition 
that  is  difficult  to  violate. 

The  lti  P-j  (s)  is  precisely  equivalent  to  the  nonlinear  w-1  (from  the 
input-output  viewpoint) ,  only  when  the  input  i s  x-j  (o£  the  output  i s  y  =  a . ) . 
Thus,  imagine  each  is  put  inside  separate,  externally  identical  black  boxes. 

The  only  measurement  you  are  allowed  is  to  inject  an  input  signal  x-j(t)  . 

It  will  be  impossible  to  say  which  black  box  contains  the  nonlinear  wJ 
and  which  one  the  linear  P-!  . 

i 

Repeat  this  operation  over  all  the  w  e  It),  for  the  same  a.  .  Then 

i  A 

repeat  it  for  all  the  a  e  A  ,  generating  a  set  of  (P-j)  =  P  .  For  example, 
if  the  set  A  had  10  elements  and  W  had  20,  then  the  set  P  would  have 
200  members.  In  practice  both  A  and  W  are  uncountable  sets  and  so  is  P  . 
The  set  P  is  equivalent  to  W  only  with  respect  to  the  set  A  .  Thus, 
for  any  pair  am  e  A  ,  wn  e  W  ,  there  exists  a  P^  e  P  which  is  equivalent 
to  wn  ,  in  the  sense  that  some  signal  must  be  applied  to  wn  to 

J. 

Implicit  in  this  is  the  assumption  that  the  plant  has  a  unique  inverse, 
excluding  for  example  hard  saturation.  In  practice,  in  such  cases,  one  can 
model  the  hard  saturation  by  a  very  small  gain.  Actually,  the  theory  can  be 
expanded  to  the  case  where  there  is  no  unique  inverse  but  there  is  a  set  of 
possible  inputs  giving  the  same  output,  providing  this  set  satisfies  certain 
reasonable  conditions. 
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its  output  is  also 
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m>  and  if  this  same  signal  is  applied  to  p”  , 

a  .  P  is  the  linear  time-invariant  equivalent  (LT IE)  ot  W,  with 
m 

respect  to  A  . 

In  Fig.  3~la,  suppose  w  (which  can  be  any  member  of  the  nonlinear 
set  W  )  is  replaced  by  P(s)  ,  any  member  at  all  of  P  .  Consider  the 
problem  of  choosing  F,  G  so  as  to  guarantee  that  the  system  output  lies 
in  the  OK  set  A,  no  matter  which  P  e  P  happens  to  be  the  plant.  This  is 
purely  a  quantitative  Iti  design  problem  for  which  the  techniques  of  Chapter  2 
may  be  applicable.  Suppose  they  are,  which  means  that  F,  G  are  found  such 
that  the  output  is  in  A  ,  no  matter  which  P  e  P  is  chosen.  Then,  under 
very  general  conditions  this  same  F,  G  compensation  pair  works  for  the 
nonlinear  set  W  .  This  means  that  no  matter  which  nonlinear  plant  w  e  W 
is  used,  the  output  is  guaranteed  to  be  in  the  OK  set  A  .  Functional 
analysis  techniques  are  used  [3,4]  to  prove  the  above,  but  the  design  execu¬ 
tion  involves  simple,  direct  frequency  response  techniques,  which  can  be 
performed  by  any  reasonably  competent  feedback  design  engineer.  A  brief 
outline  of  the  procedure  is  next  presented. 

3.3  Outline  of  Synthesis  Procedure  for  Nonlinear  Uncertain  Plants 
3.3-1  Problem  Statement 

To  simplify  the  presentation,  take  an  example  where  the  nonlinear  plant 
is  given  by  the  first  order  equation 

y  +  Ay2  =  Bx  ;  A  e  [-1,5]  ,  Be  [1,10]  (3-3) 

Thus,  the  nonlinear  plant  is  clearly  "unstable"  for  a  subset  of  plant  parameter 
values.  Suppose  the  typical  command  inputs  are  steps  from  k|  *  1  to  5 
magnitude,  and  other  inputs  of  the  form  k2te  at  with  k^  e  [1,4]  , 
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a  e  [.2,3]  •  It  is  desired  that  the  closed-loop  system  behave  like  a  1 t i 
system  with  transfer  function 

T(s)  =  -5— 1 = -  ,  E  £  [16,20]  ,  K  e  [9,11]  (3-4) 

s4  +  Ks  +  E 


3.3.2  Design  Procedure 

Take  any  command  input  r(t)  in  the  set  (kju(t)  ,  k2te  }  and  any 
acceptable  T(s)  .  Their  transform  product  is  R(s)T(s)  whose  inverse 
transform  gives  a  possible  output  y(t)  .  For  example,  if  R(s)  =  2/s  ,  then 


y (t)  =  L'’r(s)T(s)  =  2  +  2+  y(t) 

d  ”  a 

a  +  b  =  K,  ab  =  E. 


(3.5) 


Substitute  y(t)  into  Eq.  (3*3)  and  solve  for 

Bx  ( t)  =  ab(e--b-~a£-btj--*-  A  [4  +  4y(t)  +  p2(t)  ] .  (3.6) 

Next  find  P(s)  =  ■  .  Actual  ly,  it  is  easier  to  find  p^~y  ■  y|^y  wi  th 

Y(s)  =  2E/s (s2  +  Fs  +  E)  , 


BX(s) 


(5+iTC5Vby  +  A 


4(s+a+b)  s^  +  3s(a+b)  +  2  (a^+3ab+b2) 
(s+a) (s+b)  (s+2a) (s+2b) (s+a+b) 


(3-7) 


The  above  g  i  ves  a  set  (P  (s) }  ,  by  letting  a,  b,  A,  B,  kj,  etc. 
range  over  their  permissible  values.  The  process  is  repeated  by  taking 

f 

r(t)  =  l^te  ,  giving  in  the  same  manner  another  set  of  P(s)  .  The  union 
of  the  two  sets  gives  the  total  P  set.  Thereafter,  the  design  procedure 
is  entirely  the  same  as  for  an  uncertain  linear  time-invariant  P  set-- as 
briefly  described  in  Section  2. A  and  in  much  more  detail  in  [7].  The 
objective  is  to  find  F(s)  ,  G(s)  in  Fig.  3"la,  such  that  the  system  transfer 
function  satisfies  the  frequency-domain  tolerances  given  by  (3-4),  for  ail 
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?(s)  £  P  .  This  is  completely  a  linear  time  invariant  design  problem.  If 
it  is  solvable,  then  the  same  F(s)  ,  G(s)  are  guaranteed  to  work  for  the 
original  uncertain  nonlinear  plant  problem.  In  this  case  it  is  certainly 
solvable. 

3.3.3  A  Condition  for  Guaranteed  Solvability  of  the  Linear  Time  Invariant 
Equivalent  (LTIE)  Problem 

If  the  LTIE  set  P  =  (P(s)}  is  minimum-phase  then  the  LTIE  problem  is 
solvable.  Now  P(s)  =  Y(s)/X(s)  with  x(t)  =  w  '[y(t)]  in  the  notation 
of  Section  3- 1 •  Zeros  of  P(s)  are  the  zeros  of  Y(s)  (those  uncancelled 
by  X(s)  )  and  poles  (uncancelled)  of  X(s)  .  If  the  command  input  is 
minimum-phase,  then  the  OK  output  set  A  can  be  prescribed  minimum-phase, 
so  Y(s)  has  no  right  half-plane  zeros.  The  condition  then  is  that  X(s) 
has  no  right  half-plane  poles,  i.e.  the  plant  input  is  "stable." 

Consider  a  nonlinear  plant  described  by  a  differential  equation  of  the 

form  N j y  •=  N2X  with  Nj  ,  N2  nonlinear  differential  operators.  it  is 

assumed  y(t)  is  "stable"  and  so  is  N ^ y  .  For  example  all  r^y  of  the 

form  E  f. (y)  — %■  with  f .  (y)  bounded  will  give  "stable"  N.y  .  So  the 
'  dt'  '  1 

problem  is  the  stability  of  x  in  N2x  =  i/>(t)  with  ( t )  =  Njy  known  and 
"stable".  For  example,  if  N^x  =  x  -  x  (linear  in  this  case),  then  the 
plant  is  nonminimum-phase,  as  it  would  be  even  if  N j y  was  linear  time- 
invariant.  The  above  gives  a  relatively  simple  means  for  determining  whether 
the  LTIE  plant  set  is  minimum-phase,  a  concept  which  is  now  also  meaningful 
in  nonlinear  feedback  design. 
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.3.4  System  Response  to  Other  Inputs 


It  is  important  to  emphasize  that  the  desired  performance  can  be 
guaranteed  only  for  those  system  inputs  which  were  explicitly  listed  and 
considered  in  the  design.  Thus,  in  the  example  of  this  section,  what  can 
be  said  of  an  input  r(t)  =  5 t  ?  Apriori  not  very  much,  because  it  was  not 
in  the  list  for  which  the  ELT1  plant  set  was  found.  We  could  however,  check 
if  the  closed-loop  system  will  behave  like  a  T(s)  satisfying  (3-4)  to 
such  an  input,  as  follows. 

Let  Y(s)  =  R(s)T(s)  ,  with  R(s)  =  5/s2  in  this  case  and  T(s)  of 
(3-4).  Solve  (3-3)  for  x(t)  and  find  P(s)  =  Y(s)/X(s)  .  If  P(s)  e  P 
the  set  previously  obtained,  then  indeed  the  system  behaves  like  T(s)  to 
this  input  too.  However,  if  ramps  of  a  certain  size  are  indeed  typical 
inputs,  then  it  would  be  a  good  idea  to  include  them  at  the  outset  in  the  set 
of  r(t)  ,  for  which  the  design  is  explicitly  made.  In  any  significant 
problem  (such  as  the  flight  control  one  described  in  Chapter  4),  one  must 
prepare  a  computer  program  for  solving  the  nonlinear  plant  equation  backwards 
for  the  input  x(t)  ,  given  the  output  y(t)  .  A  computer  program  is  also 
prepared  f«r  finding  Y(s)  =  T(s)R(s)  .  It  is  then  a  simple  matter  to  err 
on  the  safe  side  and  let  the  set  (r(t)}  include  all  conceivable  realistic 
inputs. 

In  the  preceding,  only  system  command  inputs  r(t)  in  Figure  3“la  were 
considered.  But  disturbance  inputs  dj  and/or  d2  in  Figure  3"la  are 
similarly  treated.  One  must  assign  acceptable  output  sets  in  response  to 
the  disturbances,  in  the  same  manner  as  the  set  { r ( t ) }  .  This  gives  a  set 
of  acceptable  outputs  (y^(t)}  ,  due  to  the  disturbances.  The  corresponding 
plant  inputs  and  LTIE  plant  set  are  found  in  the  same  manner.  An  example 


in  which  both  command  and  distrubance  inputs  were  thus  considered  is  given 
in  [12].  Incidentally,  it  has  been  shown  [A]  that  the  design  is  not  sensitive 
to  the  inputs,  i.e.,  significant  deviations  from  the  r(t)  for  which  the 
design  was  made,  will  have  small  effect  on  the  T(s)  presented  by  the 
system.  The  design  can  be  executed  such  that  noninfinitesimal  changes  in 
each  r(t)  give  a  T(s)  within  the  acceptable  set.  This  has  been  verified 
in  ail  our  design  examples. 


There  are  quite  a  few  significant  extensions,  briefly  noted  here.  One 
is  to  nonlinear  closed-loop  synthesis,  i.e.  the  closed-loop  system  can  be 
designed  to  be  nonlinear  and  quantitatively  so.  For  example,  it  may  be 
desirable  to  have  the  plant  operate  close  to  its  saturated  value  nearly  all 
the  time,  in  order  to  have  fastest  possible  response,  say  to  step  functions. 
The  OK  response  sets  may  then  be  as  shown  in  F i c .  3-2.  Note  that  in  a  linear 


system  the  OK  response  set  for  a  step  input  of  magnitude  2  must  be  twice 
that  for  a  step  input  of  1.  This  is  not  so  in  Fig.  3“2.  For  any  step  input 
in  a  specified  range,  the  output  is  required  to  climb  at  the  same  (saturated) 
rate  until  it  is  close  to  the  final  commanded  value.  F  in  Figure  3_1  then 
emerges  nonlinear.  A  design  example  of  this  kind  is  described  in  detail  in 
[121.  Another  such  nonlinear  closed- loop  design  philosophy  may  be  applicable 
to  flight  control.  In  a  fast  maneuvering  situation,  the  command  inputs 
are  fast  and  large,  and  it  is  desired  to  command  acceleration  while  in  a 
tracking  mode,  the  command  inputs  are  slow,  and  it  is  desired  to  command 
pitch.  One  could  pose  typical  inputs  of  the  first  kind  and  define  the 
desired  acceptable  set  and  repeat  for  inputs  of  the  second  kind.  The  design 
technique  can  handle  this  kind  of  situation — even  with  the  plant  nonlinear 


and  uncertain,  of  course. 


The  plant  can  be  linear  or  nonlinear  time-varying  and  uncertain,  e.g. 

f  _  ry  £  \  —  f 

.  A+Be  '  .  ,  „  (H+je  )  „.  .  „  -psinvt  3  i-, 

y  +  yy  +  Qy  J  -  Kx  +  Me  (3- 


with  say  uncertainties  E  e  [1,5]  ,  G  e  [2,4]  ,  a  e  [1,2]  ,  etc.  The  closed- 
loop  system  can  be  designed  to  be  linear  time  invariant,  linear  time-varying, 
nonlinear  time  invariant  or  nonlinear  time-varying.  If  linear  time-varying 


then  F,  G  in  Figure  3*1  emerge  linear  time-varying.  If  nonlinear  time- 
varying,  F  is  nonlinear  time-varying,  G  linear  time-varying.  One  example 


of  a  highly  uncertain  nonlinear,  time-varying  plant  but  with  linear  time 
invariant  closed  loop  specifications,  is  given  in  [**].  Another  example  with 
a  linear  time-varying  uncertain  plant  and  Iti  closed-loop  specifications  is 
given  in  [3]- 

Another  extension  is  to  single  input-output  multiple-loop  systems  of 
the  kind  described  in  Chapter  2,  but  excluding  the  plant  modification  structure 
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Consider  Figure  2- 10,  with  nonlinear  replacing  P]^ . 

If  c(t)  is  known  (a  member  of  the  OK  set),  one  can  solve  the  nonlinear 
equations  backwards  to  find  C2 (t)  ,  c^(t)  ,  etc.  Then  C(s)/C2(s)  =  Pj 
the  Iti  equivalent  of  Wj  ,  C2(s)/Cj(s)  =  P2  the  equivalent  of  w 2  etc. 

A  set  of  each  is  generated  and  the  result  is  a  multiple-loop  Iti  problem 
solvable  by  the  methods  of  [19].  The  resulting  Iti  design  works  for  the 
nonlinear  uncertain  plant  set. 

Another  extension  being  currently  researched  is  for  the  case  of  non¬ 
minimum-phase  system  inputs.  Nonminimum-phase  OK  output  sets  must  then  be  postu 
lated,  leading  to  Iti  plant  equivalents  which  may  have  both  zeros  and  poles 
in  the  right  half-plane.  It  appears  that  in  the  resulting  Iti  design,  there 
can  be  zeros  of  I  +L(s)  in  the  right  half-plane,  where  L(s)  is  the  Iti 
equivalent  loop  transmission.  In  a  real  Iti  system  this  would  give  an 
unstable  system,  but  it  need  not  be  so  in  a  Iti  equivalent  system.  This 
research  has  not  as  yet  been  completed. 


3-5  Comparison  with  Other  Synthesis  Techniques 

There  are  no  other  techniques  which  cope  di rectly  wi th  the  quanti tati ve 
synthesis  problem  described  here,  i.e.  which  guarantee  outputs  in  the  time 
or  frequency  domain,  for  apriori  given  uncertainty  ranges  of  plant  parameters, 
even  For  linear  time  invariant  plants.  Optimal  control  theory  is  mostly 
concerned  with  perfectly  known  plants.  If  parameter  uncertainty  is  considered, 
it  is  usually  done  in  a  qualitative  manner.  And  system  performance  is  always 
formulated  by  a  scalar  funct iona I -- a  single  real  number  which  is  usually 
a  quadratic  function  of  the  state  and  control  variables.  In  most  control 
problems,  one  is  vitally  interested  in  transient  response,  which  is  a  function 
of  time,  not  a  scalar  functional.  Often  systems  which  are  optimal  for  the 
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quadrat i cal ly  based  functional,  have  very  poor  (large  overshoot)  in  their 
time  response. 

In  the  last  decade,  modern  control  theory  has  concentrated  on  state- 
variable  feedback,  use  of  observers,  etc.  for  realizing  desired  system 
eigenvalues.  Again,  plant  uncertainty  has  until  recently,  been  completely 
neglected.  Lately  there  has  been  a  great  deal  of  work  on  "robustness", 
mainly  with  the  problem  of  stability  for  small  enough  parameter  variations, 
or  for  disastrous  failures  o'  some  components.  This  is  a  very  welcome 
improvement  in  the  realism  t.f  modern  control  theory.  However,  it  has  very 
far  to  go  as  yet.  The  role  of  loop  bandwidth,  or  its  time-domain  equivalent 
of  speed  of  loop  response,  is  still  not  appreciated.  The  system  performance 
is  still  judged  by  a  scalar  functional.  This  laudable  attempt  at  realism 
in  control  of  uncertain  systems,  is  being  handicapped  by  the  mathematical 
formulations  and  techniques  used.  These  were  natural  and  sensible  in 
optimal  control  of  perfectly  known  plants,  which  is  basically  an  open-loop 
problem.  They  are  unnatural  and  opaque  for  quantitative  synthesis  of  highly 
uncertain  systems. 

At  the  present  time,  any  other  synthesis  technique  applied  to  the 
quantitative  problem  formulated  here,  must  proceed  by  cut  and  try.  It  is 
not  even  possible  for  any  of  them  to  declare  apriori  whether  a  specific 
quantitative  performance  specification  set  can  be  satisfied  or  not,  even  for 
lti  problems.  It  is  possible  to  do  so  with  the  techniques  presented  here, 
even  for  nonlinear  time-varying  problems.  The  next  chapter  presents  a 
detailed  application  of  the  nonlinear  technique  to  the  longitudinal  flight 
control  problem,  based  on  a  simplified  but  highly  nonlinear  plant  model. 


-46- 


Ml 


CHAPTER  4 

FLIGHT  CONTROL  DESIGN  BASED  ON 
NONLINEAR  MODEL  WITH  UNCERTAIN  PARAMETERS 


4. 1  Statement  of  Problem 

In  this  chapter,  the  breakthrough  of  Chapter  3  is  applied  to  a  signi¬ 
ficantly  nonlinear  model  of  the  short-period  longitudinal  flight  control 
problem. 


Nonlinear  plant  model  (see  Table  of  Symbols) 
c  PV0S 

u  =  qV0  +  g  cos  0  -  [^(ct)  +  CN6(a)fi] 

••  .  PV0SC  r 

6  “  q  “  ~~T\~  [Cm6(a)5  +  CMct(a)  +  W~Q  Cmq(a)q]  ’ 


(4.1) 


a  =  tan 


y 

-1  u 


V  ‘ 

vo 


(4.2,  3) 


The  control  input  is  the  elevator  angle  6  ,  and  the  output  variable  to 
be  controlled  is  [24], 


c*  =  12.4  0  +  (Ci  -  VQ0  +  69) 


(4.4) 


The  feedback  structure  used  is  shown  in  Fig.  4-1  with  c  ■  w(6)  . 

The  numbers  used  in  Equations  (4.1-4)  are  [25,26]  I  ■  207,000  kg  m  , 
m  =«  17,000  kg,  C  =  4.89  meters,  S  ■  49.2  m^.  The  C.j(a)  are  nonlinear 
functions  of  a,  see  Fig.  2-2.  Since  at  (Fig.  4-8)  ranges  in  [0,35°], 
there  is  strong  nonlinear  operation.  The  horizontal  velocity  v  was  taken 
as  VQ  fixed,  which  is  incorrect  for  some  low-velocity  cases,  but  the 
objective  here  is  to  demonstrate  the  validity  of  the  design  technique  in  a 
strongly  nonlinear  situation,  which  is  achieved  sufficiently  by  means  of 
the  nonlinear  C.j(ot)  . 
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The  bounds  on  the  acceptable  c'(t)  in  response  to  a  unit  step 
command  are  included  in  Fig.  A-7,  which  also  include  design  simulation 
results.  The  set  of  command  inputs  R  consists  of  steps  l  to  5  in 
magnitude.  Parameter  uncertainty  is  due  to  p  ranging  in  [.3  ,  1.22] 
and  Vq  in  [75.206].  Initial  conditions  are  w(0)  =  q(0)  *  q(0)  =  0, 
a(0)  =  0(0)  ,  giving  initial  values  for  6  (as  well  as  u  ,  0  ,  a)  which 
is  subtracted  out  so  that  the  change  in  6  is  used  to  find  the  1 t i 
equivalent  P  set.  The  detailed  steps  in  the  design  implementation  are 
next  presented,  with  comments  postposed  to  the  end. 

A. 2  Design  Execution 

A. 2.1  The  linear  time  invariant  (lti)  equivalent  set  P 

Let  jjc'(t)  =  C'(s)  =  T(s)R(s)  ,  where  T(s)  is  the  equivalent 
system  transfer  function  presented  by  the  closed-loop  system  of  Fig.  A-l. 
Here  R(s)  =  k/s  with  k  e  [1,5]  ,  while  it  is  required  that  T(s)  e  T, 
a  set  derived  from  the  bounds  in  Fig.  A-7.  A  simple  means  for  generating 
T,  taken  from  [25],  is  to  let  T(s)  =  a2(s + 2.9)/2.9(s2 + 2£as + a2)  ,  £ 

ranging  in  [3*7  ,  l  «5j  ,  a  in  [3.1  A  ,  7*6]  ,  giving  the  bounds  or.  |T(ju>)  j 
in  Fig.  A-3-  Such  bounds  suffice  [7]  for  mp  (and  obviously  stable) 

T(s)  .  Anv  T(s)  ,  F (s )  pair  thus  generates  an  acceptable  c  (t)  .  A 
computer  program  solved  Eqs.  (A. 1-A)  backwards  for  6(t)  and  then  checked 

JL 

the  result  by  solving  Eqs.  (A.l-A)  forwards  for  c  (t)  from  6(t)  .  The 
program  was  considered  adequate  only  when  excellent  agreement  was  obtained 
over  the  entire  range  considered.  C*(s)  was  a  priori  available  and 

A ( jto)  =  I6(t)j.  was  obtained  by  numerical  integration.  As  c  is  a 

JU3 

short-period  criterion  and  for  all  the  acceptable  cases  has  definitely 
reached  steady-state  in  A  seconds  (see  Figs.  A-7),  Eqs.  (A.l-A)  were  solved 
only  for  l  e  [G,A  seconds]  ,  and  the  constant  6(A)  was  used  for  t  >  A  . 


J  JP  *  ™ 


Loci  of  six  P(joj)  are  shown  in  Fig.  4-4,  two  of  them  (e,f)  unstable 
with  a  pair  of  right  haif-piane  poles,  which  are  zeros  of  A(s)  .  The  set 
includes  a  large  number  of  such  unstable  1 t i  P(s)  ,  which  the  design 
technique  can  easily  handle. 

Plant  templates.  At  any  oi  say  ui  =  oij  ,  the  set  {P ( joo)  >  ,  P  e  P 
consists  of  a  region  in  the  logarithmic  complex  plant  (Nichols  chart) 
denoted  as  the  (Oj-plant  template  Tp(o)j)  .  A  number  of  Tp(oj)  are  shown 
in  Figs.  4~5a  to  f.  At  very  small  U)  there  are  two  almost  constant  angle 
sub-tempi ates  360°  apart.  This  is  due  to  the  presence  of  both  stable  and 
unstable  P  e  P,  and  the  fact  that  Arg  P  near  a)  =  0  is  either  =  0 
or  nir/2  for  some  integer  n  .  As  to  increases,  the  two  groups  merge 
together  and  approach  a  single  vertical  line  at  large  w  well  beyond  the 
plant  "dynamics"  (w  =  12  is  large  enough  in  this  case,  see  Fig.  4-5f). 

Note  how  this  frequency  response  approach  is  indifferent  to  system  order. 

4.2.2  Design  of  G(s)  ,  F(s) 

Given  the  set  P  =  (P)  ,  the  problem  is  to  find  F(s)  ,G(s)  in 
Fig.  4-1  such  that  the  system  transfer  function  T(joj)  *  FGP/O+GP)  e  T, 
for  all  P  e  P  .  One  may  program  the  computer  to  find  the  (unique)  bounds 
B  (to)  on  G(jw)  ,  so  that  as  P  ranges  over  P,  A  1  n  |  T ( jco)  j  = 

GP 

A  1  n  |  |  ^^(u))  -  Aj(td)]  db  of  Fig.  4-3.  Alternatively,  this  may 

be  done  by  hand,  giving  useful  insight  as  explained  in  Section  2.4.3, 

Figure  2-4. 

Fig.  4-6  shows  the  bounds  so  obtained  on  G(ju>)  ,  and  the  G(jw) 

chosen  to  satisfy  these  bounds.  Let  eg  be  the  excess  of  poles  over  zeros 

eG 

assigned  to  G(s)  ,  so  that  as  s  -*■  »  ,  G  -*•  kg/s  .  It  is  reasonable  to 
define  the  optimum  G  as  that  which  satisfies  its  bounds  with  minimum  kg  . 
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It  has  been  shown  [11]  that  G _ .  lies  on  B(u)  at  all  a>  and  that  G 

opt  opt 

exists  and  is  unique.  The  design  of  a  practical  G(s)  to  satisfy  the 
bounds  is  somewhat  of  an  art  [7].  For  a  given  skill  in  the  art,  the  greater 
the  number  of  poles  and  zeros  of  G  ,  the  closer  one  can  get  to  the  optimum, 
so  there  is  trade-off  between  complexity  and  bandwidth.  Here,  we  chose 
simply  by  cut  and  try  G(s)  =  (1  +  .2s)/s(l  +  .033s2)  (1  +  .002s  +  10  **s2) 
with  very  modest  bandwidth,  Figs.  4-3,  6.  A  much  simpler  G(s)  could  have 
been  chosen  with  larger  bandwidth.  The  designer  must  make  his  own  trade¬ 
off.  Ref.  7  offers  some  advice  on  the  shaping  of  a  function  to  satisfy 
a  set  (B(u))]  in  the  Nichols  chart. 

Design  of  F (s) .  G (s )  only  guarantees  that  A | T ( joi)  |  £  A„  (w)  -  Aj  (w)  of 
Fig-  4-3,  e.g.,  at  oo  =  10,  the  actual  change  in  | L ( j  10)/ ( 1  +  L ( j  10) )  |  is 
from  -7  db  to  4  db  ,  while  from  Fig.  4-3,  the  permitted  change  in 
I T ( j  1 0)  ]  *  |FL/(1+L)|  is  from  -15  to  2.8  db  .  Hence,  any  value  of 
I F  ( j 1 0) |  e  [-8,  -1.2  db]  is  acceptable.  In  this  way,  upper  and  lower 
bounds  on  [  F  ( ju>)  ]  are  obtained  and  F(s)  is  chosen  to  satisfy  them, 
which  is  also  somewhat  of  an  art.  In  this  example,  a  satisfactory 
F (s)  =  (1  +  .33s) (1  +  .05s)/(l  +  .25s) (1  +  .2s) (1  +  .0125s)2  ,  see  Fig.  4-3. 

4.2.3  Design  Simulation 

The  nonlinear  system  was  simulated  and  its  response  found  for  several 
hundred  command  inputs  and  gust  disturbances.  Typical  (68)  responses 

J. 

to  c  step  commands  are  shown  in  Fig.  4-7  a  to  h  for  various  combinations 
of  Vq  ,  p  and  step  (k)  values.  The  transient  response  of  a(t)  ,  0(t) 
etc.,  depend,  of  course,  on  the  values  of  k  ,  p  ,  VQ  .  Two  sets  of  these 
are  shown  in  Figs.  4-8  a,  b,  with  Fig.  8a  depicting  very  large  a(t) 
excursion,  for  which  the  C.j(a)  in  Fig.  4-2  are  in  strongly  nonlinear 
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Figures  4-7  c,  d.  Simulation  results.  Response  to  step 
command  of  c  ,  magnitude  K. 
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Figures  k~7  q.  h.  Simulation  results.  Response  to  step 
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command  of  c  ,  magnitude  K. 
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ranges,  These  are  the  inputs  for  which  the  system  was  designed,  and  for 
which  guarantees  can  be  made.  Out  of  the  several  hundred  cases  simulated 
in  a  very  few  cases  (Fig.  4-7  h) ,  there  was  very  slight  excursion  out  of 
the  bounds.  It  is  possible  to  include  in  the  design  other  inputs  and  gust 
disturbances,  with  specified  response  tolerances  -  and  then  guarantees 
can  be  made  for  these  3s  well.  The  response  to  other  inputs  is  nevertheless 
found  here  to  be  also  quite  satisfactory.  This  is  typical  of  the  design, 
i.e.,  the  system  is  not  very  sharply  tuned  to  the  class  of  inputs  used  in 
the  design  execution.  There  is  reasonable  response  continuity  to  other 
inputs. 

Some  responses  to  very  large  c  step  commands  causing  hard 

saturation  are  shown  in  Fig.  4-9  a,  b.  Response  to  Gust  Disturbances. 

The  gust  input  was  modelled  by  replacing  a  in  (4.3)  by  a  =  tan  *  r- +  a 

*0  9  s 

Two  kinds  of  a  _  were  used.  In  one  a  .  is  a  half-sine  wave  of 
gust  gust 

amplitude  20/Vq  radians  and  half-period  THALF  z  [.2,2]sec  .  Some 
results  are  shown  in  Figs.  4-10  a  to  c.  In  (a)  the  gust  begins  precisely 
at  the  instant  of  application  of  simultaneous  c  step  commands.  The 
second  kind  in  Fig.  4-11  a-d,  is  stochastic  guassian  with  power  spectrum 
k/(l+w  )V0  and  (agUSt)rms  =  6/V ^  radians.  Examples  of  responses  to  a 

A 

single  square  wave  c  command  with  equal  positive  and  negative  values  k 
and  total  duration  2  THALF  are  shown  in  Fig.  4-12  a  to  c. 
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Figure  4-10  b.  Responses  to  gust. 
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4.3  Discussion 

A  second-order  model  was  used  for  T(s)  with  excess  of  poles  over 
zeros  ey  ■  1  ,  in  Step  1.  This  appears  to  be  incompatible  with  G,  F 
fifth  order  and  excesses  eg  ■  4  ,  ep  *  2  ,  and  P  (no  analytical  expression 
but  linearized  (4. 1-4)  give  second-order  P)  with  ep  »  0  .  Strict  design 
execution  appears  to  require  a  T(s)  of  complexity  compatible  with 
T  =  F  GP/O+GP)  .  Note,  however,  that  in  Step  1,  the  model  of  T(s)  is 
used  only  to  generate  a  set  which  covers  the  range  of  a  priori  specified 
acceptable  outputs.  Any  T(s)  model  which  achieves  this  is  clearly 
satisfactory,  and  the  simpler  the  better.  The  designer  can  later  choose 
the  complexity  of  G(s)  ,  F(s)  ,  with  no  regard  for  that  of  the  T(s) 
model  used  in  Step  1. 

The  class  of  applicable  non 1  inear i ti es  has  been  defined  implicitly  in 
Ref.  4,  but  one  very  large  class  can  be  defined  explicitly.  Let 
Djy(t)  =  D^xCt)  with  D|  ,  D2  operators  which  may  be  nonlinear,  uncertain 
and  time-varying,  e.g. 

— ot  t 

D,y  =  M(y)(y)‘5sgn  y  +  -yyyt —  y2  I  y  I P  +  Hy2 

Ee  [1,5,  Fe  [.5,4],  A  e  [-3,6],  ae  [-5,1.5] 

B  e  [-3,2]  ,  H  e  [-4,1]  ,  n  0  [.5,2]  ,  M  r  [1,5]  • 

The  range  of  M  must  be  of  the  same  sign.  A1 1  OK  y  and  D1 y  must  be 
bounded  for  all  t  in  [0,°°]  and  D^y  must  exist.  Hence  y  must  be 
twice  differentiable  except,  at  most,  at  a  countable  number  of  points. 

Then  D j y  =  ip ( t )  is  known  and  there  must  exist  a  unique  solution  for  x 
in  D2X  =  ip ( t )  .  The  solution  must  be  bounded  for  all  t  £  [0,°°]  .  Thus, 

0^*  *  y(t;  must  be  "bounded- input ,  bounded-output"  stable.  However, 

^,/  ti  -  ;•(  t)  may  be  "unstable"  in  that  a  bounded  v(t)  is  allowed  to 


-71- 


result  in  unbounded  z(t)  .  It  is  only  necessary  that  bounded  z(t)  gives 
bounded  v(t)  . 

It  is  possible  that  a  simple  linearization  might  do  just  as  well  in 
practice,  but  this  is  a  matter  of  chance,  whereas  this  design  technique  is 
guaranteed  to  work  if  the  constraints  are  satisfied. 

The  constraint  on  W  that  P  is  mp,  is  required  because  only  then 
can  one  guarantee  that  any  specifications  no  matter  how  narrow  (A2  “  A. 
arbitrarily  small  but  nonzero  in  Fig.  4-3),  may  be  satisfied  for  arbitrarily 
large  but  bounded  parameter  uncertainty  (but  some  parameters  must  not  change 
sign  [23]).  No  such  guarantee  can  be  made  for  nmp  P,  but  the  problem 
is  still  solvable  if  the  specifications  are  not  too  narrow  and  P  is  not 
too  large  [II]  a  set. 
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CHAPTER  5 

A  SECOND  BREAKTHROUGH  IN 
QUANTITATIVE  FEEDBACK  SYNTHESIS 

5.1  The  Muliple  Input-Output  (mio)  Problem 

A  very  large  number  of  control  problems  involve  multiple  inputs  and 

outputs,  with  each  control  input  affecting  all  the  output  variables  to  some 

extent.  This  is  certainly  so  in  many  flight  control  modes.  In  the  n  x  n 

2 

case,  there  is  a  matrix  of  n  system  response  functions  to  realize.  In 
realistic  control,  the  n  x  n  plant  matrix  has  significantly  uncertain 

parameters.  The  rea I  control  problem  is  a  quantitative  one:  to  satisfy 

2  2  2 
specified  2n  tolerance  sets  on  the  2n  response  functions  (n  to 

2 

command  inputs  and  n  to  disturbance  inputs),  despite  the  specified  uncer¬ 
tainties  in  the  parameters  of  the  n  x  n  plant  relations. 

There  is  a  vast  literature  on  this  mio  (often  denoted  as  "multivariable") 
problem.  However,  nearly  all  of  it  deals  with  assumed  perfectly  known 
plants,  so  the  problem  treated  is  securing  the  desired  closed  loop  response, 
under  the  constraint  of  a  feedback  structure  around  the  perfectly  known 
plant.  Just  as  in  the  single  input-output  case,  the  quantitative  aspects -- 
the  extent  of  the  plant  uncertainties  and  the  performance  tolerances  (i.e., 
the  very  reasons  for  using  feedback)  --do  not  enter  as  design  parameters. 
Stabilization  of  the  resulting  highly  mul t Iple- loop  system  is  considered 
a  major  undertaking.  But  note  that  this  is  being  done  for  a  perfectly 

known  plant  matrix.  Consider  how  much  more  difficult  is  this  stabilization 

2 

problem  when  the  parameters  of  the  n  plant  functions  have  considerable 
uncertai nty -- not  10%  or  20%,  but  hundreds  or  thousands  per  cent.  Furthermore, 
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such  stabilization  is  only  one  part  of  the  problem  in  genuine  quantitative 

synthesis.  Stabilization  is  insufficient,  for  one  must  also  guarantee  that 

each  of  the  2n  closed  loop  input-output  functions  satisfy  specified 

tolerances,  over  the  range  of  plant  parameter  uncertainty. 

Thus,  the  quantitative  Iti  mio  problem  with  large  uncertainty  has  long 

appeared  to  be  intractable,  except  by  cut  and  try.  There  has  recently 

appeared,  however,  an  exact  synthesis  technique  for  this  problem  [23; 

Appendix  I  of  this  report].  A  highly  attractive  feature  of  this  technique 

is  that  for  the  n  x  n  case,  it  consists  of  designing  n  distinct, 

2 

separate,  non  interacting  loop  transmission  functions  and  n  prefilter 
functions  of  the  F  type  in  Figure  3-la.  Thus,  the  n  x  n  problem  is 
broken  up  into  a  number  of  separate  single-loop  problems.  The  interaction 
between  them  is  that  of  specifications.  The  tolerances  on  the  output  of  one 
of  the  single-loop  problems,  appear  as  "disturbances"  in  a  different  one. 

It  may  be  worthwhile  to  deliberately  tighten  the  sped f icat ions  on  some  of 
the  outputs,  and  thereby  decrease  some  loop  bandwidths  needed,  without 
increasing  any  others.  Sometimes  there  is  trade-off,  in  that  such  tightening 
decreases  some  but  increases  others.  However,  these  do  not  affect  the  basic 
isolation  and  noninteraction  of  the  separate  single-loop  designs.  They 
affect  the  optimality  of  the  designs  in  terms  of  loop-bandwi dths  needed  to 
do  the  job. 

The  above  is  a  very  desirable  property  of  the  design  technique.  The 

stabilization  of  a  feedback  system  containing  a  highly  interacting,  highly 

uncertain  n  x  n  Iti  plant  becomes  one  of  stabilizing  n  distinct,  separate, 

non  interacting  single  loops.  But  even  more  than  that  —  the  quantitative 

synthesis  problem  is  solved.  The  apriori  assigned  performance  specifications 
2 

on  the  2n  response  functions  are  achieved  over  the  entire  range  of 
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parameter  uncertainty.  There  is  no  need  for  design  iteration  after  the 
separate  single  input-output  systems  have  been  designed.  This  is  for 
minimum-phase  plants.  In  the  nonmlnimum-phase  case,  the  performance  tolerances 
must  be  compatible  with  the  nonminimum-phase  character  of  the  plant,  just  as 
in  single  input-output  systems.  The  research  work  leading  to  the  above 
results  was  performed  under  AFOSR  sponsorship  [23],  but  because  of  its 
importance  in  flight  control,  is  reproduced  here  as  Appendix  I.  Since  then 
we  have  gained  some  experience  with  the  technique  and  learned  how  to  achieve 
the  trade-offs  mentioned  above. 

5-2  Extension  to  Nonlinear  Uncertain  Multiple  Input-Output  Plants 

The  nonlinear  design  philosoDhv  of  ChaDter  k  can  be  combined  with  that 
of  Section  5-1.  giving  a  quantitative  synthesis  technique  for  highly  uncertain 
nonlinear  n  x  n  multiple  input-output  plants.  The  essential  features  of 
both  synthesis  methods  are  retained.  Because  of  uncertainty,  there  is  not 
one  set  of  n  x  n  nonlinear  plant  relations,  but  an  infinite  class  of  such 
sets.  This  class  W  is  replaced  by  an  infinite  class  P  of  n  x  n  matrices 
of  linear  time  invariant  relations.  P  is  equivalent  to  W  with  respect  to 
the  sets  of  acceptable  outputs.  To  do  this  the  nonlinear  relations  must  be 
solved  backwards  to  find  the  plant  inputs  required  to  give  the  acceptable 
outputs.  The  computer  is  an  essential  tool.  Once  P  is  available,  the 
problem  becomes  a  quantitative  linear  time- invariant  multiple  input-output 
one,  for  which  the  technique  of  Appendix  1  may  be  used.  Section  5  of 
Appendix  1  presents  a  somewhat  more  detailed  outline  of  the  procedure. 
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5.3  Application  of  Nonlinear  and  Multivariable  Synthesis  Techniques  to 
Flight  Control 

In  many  of  the  Important  operating  modes,  the  flight  control  problem  is 
primarily  one  of  regulation  and  control  despite  large  uncertainty.  It  is 
one  of  achieving  desired  trajectories  as  a  function  of  time.  It  is  not  one 
of  minimizing  some  quadratic  functional  of  states  and  control  inputs.  In 
fact,  in  certain  situations  one  wants  to  drive  the  control  variables  as  hard 
as  possible.  The  justification  nevertheless  for  use  of  optimal  control 
despite  this  and  despite  its  ignoring  parameter  uncertainty  and  performance 
specifications,  is  only  that  it  provides  a  trial  solution.  But  it  is  clearly 
not  a  synthesis  tool  developed  for  the  flight  control  problem.  On  the  other 
hand,  the  quantitative  synthesis  techniques  outlined  in  Chapters  2,  3,  5 
and  illustrated  in  Chapter  4  for  a  detailed  nonlinear  problem,  are  very  close, 
almost  tailored,  to  the  flight  control  problem.  For  the  first  time  in 
feedback  control  history,  there  are  available  precise  rigorous  synthesis 
techniques  to  cope  simultaneously  with  all  the  following  important  factors 
in  the  flight  control  problem: 

1.  time  and  frequence  domain  tolerances  on  the  output  variables; 

2.  large  parameter  uncertainties; 

3-  highly  nonlinear  plant  relations; 

4.  strongly  interacting  multiple  input-output  plants; 

5-  tolerances  which  are  functions  of  the  command  inputs,  i.e.,  a  desired 
nonlinear  closed-loop  system. 
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APPENDIX 

QUANTITATIVE  SYNTHESIS  OF  UNCERTAIN  MULTIPLE 
INPUT-OUTPUT  FEEDBACK  SYSTEM1" 

★ 

Isaac  Horowitz 

ABSTRACT 

There  is  given  an  n  input,  n  output  plant  with  a  specified  range  of 

.  O 

parameter  uncertainty  and  specified  tolerances  on  the  n  system  response  to 

2 

command  functions  and  the  n  response  to  disturbance  functions.  It  is  shown 
how  Schauder's  fixed  point  theorem  may  be  used  to  generate  a  variety  of 
synthesis  techniques,  for  a  large  class  of  such  plants.  The  design  guarantees 
the  specifications  are  satisfied  over  the  range  of  parameter  uncertainty.  An 
attractive  property  is  that  design  execution  is  that  of  successive  single¬ 
loop  designs,  with  no  interaction  between  them  and  no  iteration  necessary. 
Stability  over  the  range  of  parameter  uncertainty  is  automatically  included. 

By  an  additional  use  of  Schauder's  theorem,  these  same  synthesis 
techniques  can  be  rigorously  used  for  quantitative  design  in  the  same  sense 
as  above,  for  nxn  uncertain  nonlinear  plants,  even  nonlinear  time-varying 
plants,  in  response  to  a  finite  number  of  inputs. 
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Note:  The  Appendices  to  this  paper  are  omitted  here.  They  are  available 
in  [233. 
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QUANTITATIVE  SYNTHESIS  OF  UNCERTAIN  MULTIPLE 
INPUT-OUTPUT  FEEDBACK  SYSTEMS 

I.  INTRODUCTION 

There  is  great  interest  in  multiple  input-output  (mio)  feedback  systems, 
for  obvious  reasons.  A  great  deal  of  significant  work  (too  numerous  to  list 
but  [1-10]  are  representative  and  include  bibliographies)  has  been  done, 
primarily  in  the  realization  and  properties  of  the  closed-loop  input-output 
relations,  under  the  constraint  of  a  feedback  structure  around  the  known, 
fixed  mio  "plant."  There  has  been  notable  work  done  with  uncertain  inputs, 
but  again  only  with  fixed,  known  plants.  Of  course,  plant  uncertainty  is 
always  implicit,  if  only  because  of  the  usual  approximations  required  to 
obtain  a  linear  time-invariant  model. 

In  any  case,  there  does  not  exist  as  yet  any  "quantitative  synthesis" 
technique  for  the  mio  problem  with  significant  plant  uncertainty,  even  for  the 
linear  time-invariant  case.  By  "quantitative  synthesis"  is  meant  that  there 
are  given  quantitative  bounds  on  the  plant  uncertainty,  and  quantitative 
tolerances  on  the  acceptable  closed-loop  system  response.  The  objective  is 
to  find  compensation  functions  which  guarantee  that  the  performance  tolerances 
are  satisfied  over  the  range  of  the  plant  uncertainty.  In  "quantitative 
design,"  one  guarantees  that  the  amount  of  feedback  designed  into  the  system 
is  such  as  to  obtain  the  desired  tolerances,  over  the  given  uncertainty  range. 
In  other  designs,  the  amount  of  feedback  may  be  more  or  less  than  necessary-- 
it  is  a  matter  of  chance.  The  practical  experienced  designer  may  find  the 
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latter  approach  sufficient.  However,  a  scientific  theory  of  feedback  should 
certainly  include  quantitative  design  techniques. 


In  this  paper  it  is  shown  how  Schauder's  fixed  point  theorem  can  be 
used  to  generate  a  variety  of  precise  quantitative  mio  synthesis  techniques 
suitable  for  various  problem  classes.  An  outstanding  feature  of  each  synthesis 
procedure  is  that  it  consists  of  a  succession  of  direct  (no  iterations 
necessary)  single-loop  design  steps.  Furthermore,  by  a  second  use  of 
Schauder's  theorem,  the  techniques  are  rigorously  applicable  to  quantitative 


synthesis  of  nonlinear  uncertain  mio  feedback  systems.  This  paper  concentrates 
on  existence  proofs  but  a  2  x  2  example  is  included. 

1.1  Preliminary  Statement  of  a  Linear  Time  Invariant  MIO  Problem 

In  Fig.  1,  P  =  [p . .  ( s ) ]  is  a  n  x  n  matrix  of  the  plant  transfer  functions 

•J  ,  * 

in  the  form  of  rational  functions,  each  with  an  excess  e. .  >  0  of  poles  over 

*  J 

zeros,  and  with  a  bounded  number  of  poles.  The  p..(s)  are  functions  of  q 

*  J 

physical  parameters,  with  m  an  ordered  real  q-tuple  sample  of  their  values. 

M  =  {m}  is  the  class  of  all  possible  parameter  combinations.  The  elements  of 
the  n  x  n  Iti  compensation  rational  transfer  function  matrices  F  =  [f-^(s)3, 

*  o 

G  =  (g--(s)j  are  to  be  chosen  practical  (each  with  an  excess  of  poles  over 
i  J 

zero).  They  must  ensure  that  in  response  to  command  inputs  the  closed-loop 
transfer  function  matrix  T  =  [tuy(s)]  (of  c  =  Tr)  in  Fig.  1  where  c,  -r  are  the 
n  x  1  matrices  (vectors)  of  system  outputs  and  inputs,  respectively,  satisfy 
conditions  of  the  form 

0  <  A  (oj)  s  |t  (ju>)|  s  B  (u>),VmeM 


0) 


sk  _  • 


If  the  tuy(s)  have  no  poles  or  zeros  in  the  right  half-plane  (are  stable  and 
minimum-phase),  then  tuy(s)  is  completely  determined  by  |tuv(jw)|,  so  (1) 
suffices  (Bode  1945).  It  has  been  shown  (Horowitz  1976)  that  time-domain 
tolerances  of  the  form 

“i<‘>  ^ 

dt 

v  =  0,  1,  . . . ,  n-j  any  finite  number,  can  be  satisfied  by  means  of  tolerances 
like  (1)  on  |c(ju))j,  where  c(s)  =  Jc(t).  The  writer  finds  it  much  more 
convenient  to  develop  the  synthesis  theory  in  the  frequency  domain,  and  the 
above  proves  its  sufficiency  for  time-domain  synthesis. 

This  presentation  concentrates  on  the  command  response  problem,  but  the 
same  ideas  can  be  used  to  handle  the  quantitative  disturbance  response  problem 
under  plant  uncertainty,  as  will  be  shown  in  Sec.  6.  The  constraints  on  the 
plant  and  the  specifications  are  introduced  as  needed,  in  order  to  clarify  the 
reasons  for  their  need. 


2.  DERIVATION  OF  SYNTHESIS  TECHNIQUE 
? 

In  Fig.  1,  there  are  available  n  loop  transfer  functions  in  L  = 

U, As))  =  PG,  and  n  f..  in  F  for  satisfying  the  tolerances  (1)  on  the  n 

1  J  I  J 

t. ..  But  In  the  expansion  of  T  =  [ t . - ( s ) ]  =  (I  +  l)_1LF,  each  t  .  (s,m) 

*  J  1 J  3D 

(meft)  is  a  function  of  all  the  £..(s,m)  each  uncertain,  resulting  in  very 

*  J 

complicated  expressions  for  t  ^  and  making  direct  quantitative  synthesis 
seemingly  impossible--at  least  so  far  unsuccessful.  The  objective  here  is 
to  convert  each  ^(s.m)  design  problem  into  an  equivalent  single-loop  problem 
with  uncertainty.  This  is  done  for  each  t  h,  by  lumping  all  the  other  inter- 
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acting  t. .  variables  into  an  'equivalent  disturbance* ,  as  follows. 

•  vJ 

In  Fig.  1,  c  =  PG(Fr  -  c),  so 

(P"1  +  G)c  =  GFr.  (2) 

Hence,  the  following  restriction  on  P: 

(PI):  a(s)  -  determinant  P(s)  |  O.VmeM. 


Let  ry  f  0  and  r.  =  0,  i  i1  v,  so  the  resulting  Cj(s)  =  t  j  y  ( s )  ry .  Let 

P"1  -  Cpij(s)].  (3) 

The  uth  element  of  (2)  is  then 


rv(s)  ^  ^Pui  +  9ui^tiv  |  guifiv‘ 


To  simplify  the  presentation,  we  take  gui  =  0  for  u  f  i  (although  in  practice 
it  may  be  useful  not  to  do  so).  Then  letting  ry(s)  =  1,  the  last  equation  can 
be  written  as 


uv 


uv 


JL  a  f  uv 

P  ^uuuv  P 
UU  UU  £ 

g 

1  + 

UU 


X  -  T  ,  d 

uv  duv  uv 


=  J  P  .t. 
A.  ui  i\ 


(4a) 

(4b) 


i7u 


This  corresponds  precisely  to  the  single-loop  problem  of  Fig.  2,  with 
^uve  =  l/V  course,  the  t^  in  duy  of  (4b)  are  not  known  but  the  boi 

(1)  on  |t.u|  are  known,  generating  a  set  Duv  =  { diiw } .  We  define  the  extreme  d 


1  v 1 


UV' 


lduvel  ■  «“^u  lPuiHBivl,  Bjvof  (,) 


UV 

(5) 
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Suppose  we  can  find  9UU(S)  and  fuy(s),  such  that  in  the  notation  of  (4,5) 

0  <  Ituv!  4  lTduv!!duvele  CAuv*  BUv]^meM  (6) 

Then  the  magnitude  of  the  right  side  of  (4a)  e[Auy,  Buy]  for  all  meM  and  for 
all  possible  combinations  of  t^y  (i  t  u)  which  satisfy  (1).  Suppose  this  is  so 

\/  u.v  combinations,  and  the  other  Schauder  conditions  of  Sec.  2.1  are 
satisfied.  Then  Schauder' s  fixed  point  theorem  can  be  used  to  prove  that 

O 

these  same  n  gU(J  and  n  f.y  are  a  solution  to  the  synthesis  problem  (1). 

2. 1  Application  of  Schauder' s  Fixed  Point  Theorem 

This  theorem  states  that  a  continuous  mapping  of  a  convex,  compact  set  of 

a  Banach  space  into  itself,  has  a  fixed  point  (Kantorovich  and  Akilov  1964). 

2 

We  define  the  Banach  space  to  be  the  n  C[0,°°]  product  space  denoted  here  by 
C(n^),  with  norm  =  E  individual  sup  norms.  C[0,°°]  is  the  Banach  space  of  real 
continuous  ^unctions  f(w),  we[0,°°J  with  ||f||  =  sup  |f(ui)|.  The  convex  compact 
set  in  each  of  the  n  C[0,°°]  is  taken  as  the  acceptable  set  of  |tuy(jw)|  satis¬ 
fying  (1),  denoted  by  {he(w)}  =  Huy.  Additional  constraints  have  to  be  assigned 
to  the  hg  (w)  in  order  that  each  Huy  set  is  compact  and  convex  in  C[0,<*>].  These 

constraints  have  been  justified  in  detail  in  (Horowitz  1975)  and  are  therefore 

2 

only  summarized  here.  If  each  set  is  convex  and  compact  in  C[0,°°],  their  n 

2  7 

product  set  denoted  by  H(n  )  is  convex  and  compact  in  C(n  ). 

Constraints  on  Huv  =  { h (oo) }  uv 

1.3  continuous  functions  Auv(w)>  B  (u>)  with  properties  of  (1)  as 
bounds  on  h(w) 

2.  h’(w)  is  uniformly  bounded:  3  3  |h'(w)|  <K,  V  h,u> 
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3.  h(oj)-K)  as  in  the  form  k/uie,  e  a  fixed  finite  number  >3  to  allow 

at  least  one  excess  of  pole  over  zeros  for  the  elements  of  F,G,P  in  Fig.  1 
These  constraints  guarantee  (Horowitz  1975)  that  h(aj)  can  be  taken  as  the  magni- 

A 

tude  of  a  function  h(s)s_ju  which  has  no  zeros  or  poles  in  the  interior  of 

A 

the  right  half-plane  or  on  the  jto  axis.  Arg  h(ju>)  is  obtained  from  h(w) 
by  anyone  of  a  number  of  Bode  integrals  (Bode  1945). 

An  element  of  H(n  )  consists  of  n  positive  functions  on  h;k(aj). 

Using  any  appropriate  Bode  integral,  find  the  associated  phase  function  denoted 
here  by  arg[h.-k(u>)] ,  giving  the  minimum-phase  stable  function  hik(s), 

A 

h^tjw)  =  hik(w)  +  J  argCh^w)].  For  future  use,  denote  this  sequence  of 
operations  whereby  h(w)  is  transformed  into  h(jio),  as  the  "Bode  transformation" 
B(h(w)).  Define  $  on  H(n2)  by 


*  ('I'l]  »4,]2*  *  *  *  *^*nn  ^  ^  ^uv^ll’^12*  ■*** 


9uufuv  ".^u  PuiB^hiv^w^ 


P  [1  +  IMSL) 
uu(  Puu 


(7) 


using  for  Pul- ,  Pu(j  any  specific  fixed  meM.  (Note  the  similarity  of  (7)  to(4a,b)). 
In  Appendix  2,  it  is  shown  that  guU,  fuv  can  be  found  such  that  4>  maps 

p 

H(!n  )  into  itself.  It  is  also  necessary  to  prove  4>  is  continuous,  as  follows. 


0  is  a  continuous  mapping 


? 

4>  is  continuous  if  each  of  its  n  components  is  continuous.  The  first  step 
in  each  mapping  is  B(h.v(u>))  =  (jo>).  In  (Horowitz  1975,  Sec.  Ill)  it  is  proven 
that  the  step  h^y(o))-*-  arg  v (aj)^e^ v (uj)  is  continuous  in  the  C{0,°°)  norm.  Hence, 
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the  mappings  hiv(toh  hi v(u>)  cos0.v(io)“^iv(a)) ,  hiv(u>)-*-  h1v(w)  sin  eiv(u) 
=Xiv(w)  are  continuous.  The  denominator  of  (7)  is  a  constant  on  H(n2),  and 
so  are  guu  f  and  the  Pui  in  the  numerator.  Thus,  the  numerator  has  the  form 


Num.  =  |Ka  +  jKb  -  Z(C.  +  jD.)tf.(<o)  +  jX.(u))),j  =  A, 

all  other  terms  real  and  only  thej*^,  X^  mappings  on  H(n^).  Infintesimal  changes 

in^li ,  X.  clearly  result  in  similar  change  in  Num.,  so  Num.  is  continuous  on 
2 

H(n  )  and  so  is  each  ^uv  of  (7)  and  hence  <t>.  The  conditions  in  Schauder's  theorem 
are  satisfied,  so  4>  has  a  fixed  point. 


This  means  3  a  set  of  h.  .  (w)  denoted  by  h..(w),  3 

'  J  *  J 


V“>  ■ 


^uu^uv  ^ui^iv^^ 


P  (1  + 

^  PJ 


(8) 


*"k  4: 

u , v  =  1,  ....  n,  where  h.v(jw)  =  B(hiv(u)). 


We  would  now  like  to  deduce  from  (8),  that 


B(huvM) 


g  f 
auu  uv 


- 1 
ita 


PUihUv<-’”) 


Puu<>  *  F 


uu> 


uu 


(9  ?) 


For,  if  (9)  is  true,  then  by  letting  huv(jw)  =  tuv(jio),  we  have  recovered  (4) 

o  *  + 

and  the  n  huv(ju))  are  a  solution  to  the  mio  problem  for  that  specific  meM. 
The  solution  is  unique  if  every  building  block  in  the  mio  system  has  a  unique 
output  for  any  given  input,  which  is  a  very  reasonable  condition.  This  makes 
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It  unnecessary  to  prove  that  there  are  no  transitions  from  (8)  to  an 
expression  similar  to  (9)  but  with  right  half  plane  poles  and/or  zeros.  Since 
m  is  any  element  of  M,  this  Is  true  for  all  meM  (of  course  with  a  different 
set  of  huy  for  each  m). 

The  step  from  (8)  to  (9)  is  a  crucial  one  and  must  be  justified  with  great 
care.  Given  an  analytic  function  4>(s),  there  is  an  infinitude  of  *(s)  such 
that  U(ja))|  =  |*(j»)|,  u  e[0,»],  e.g. 

(1  -  t,s)  (1  +  t?s) 

*(s)  =  *(s)  (1  +  t-j s)  (1  -  t2s) 

But  <f>(s)  /  Tp ( s )  even  though  | d> ( ja>)  [  5  |^(jto)|.  But  suppose  we  know  from  other 
sources  that  <j>^(s)  has  no  right  half  plane  zeros  or  poles,  then  given 

|+1(ju)|  =  M(u>)  a  magnitude  function  which  is  Bode  transformable,  we  can 

conclude  that  =  B(M(oj))  =  M(ju>).  Hence,  to  justify  (9)  we  must  prove 

that  the  expression  inside  the  vertical  bars  in  (8)  has  no  right  half-plane 

zeros  or  poles.  The  pole  part  is  easy,  because  1  +  9UU/PUU  is  obviously 

designed  to  have  no  right  half-plane  zeros;  certainly  guu»  fuy  won't  be 

assigned  any  such  poles;  v ( s )  doesn't  have  any  by  definition,  and  Pu^  is  not 

allowed  any  such  poles--see  Sec.  3.1.  To  prove  the  zero  part,  note  that  from 

(6)  and  Rouche's  theorem,  the  number  of  zeros  of  the  right  side  of  (9)  in  the 

right  half-plane,  equals  such  number  of 

g  f 
3uu  uv 


P  (1  +  —) 
uuu  P....; 
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which  is  easily  made  zero  in  the  single-loop  synthesis  steps  (if  Puu  has  no 
right  half-plane  poles,  a  condition  necessary  for  other  reasons— see  Sec.  3.1). 
Thus,  the  expression  inside  the  bars  in  (8)  has  no  right  half-plane  poles  or 
zeros,  justifying  (9).  This  is  a  very  valuable  result.  The  problem  of 
stabilizing  a  highly  uncertain  n  x  n  mio  system  is  automatically  disposed  of  in 
the  synthesis  procedure,  which  is  furthermore  one  of  designing  n  single-loop 
transmission  functions. 

It  is  worth  noting  that  even  if  the  above  proof  were  not  available,  it 
would  not  be  disastrous  for  this  synthesis  theory.  It  would  only  be  necessary 
to  guarantee  that  at  one  meM,  the  system  is  stable  and  minimum-phase.  For 
then,  this  would  be  soVmeM,  because  by  the  continuity  of  the  poles  (and  zeros) 
with  respect  to  the  parameters,  the  right  side  of  (8)  would  have  to  be  infinite 
(zero)  at  some  o>,  in  order  that  for  some  meM  the  system  should  be  unstable 
(or  have  a  right  half-plane  zero).  However,  the  synthesis  procedure  by 
definition  precludes  this.  And  it  is  a  relatively  easy  matter  to  guarantee 
the  desired  conditions  at  one  meM. 

3.  CONSTRAINTS  ON  MIO  PLANT 

The  above  results  hinge  on  our  ability  (a)  to  find  gyu  and  f  to  satisfy 
(6)^,  all  u,v  pairs  and  all  mqft  (b)  that  each  equivalent  single-loop  design 
is  stable  and  minimum-phase  V  mef^.  These  lead  to  constraints  on  the  mio  plant, 
obtained  by  applying  single-loop  design  theory  to  achieve  (a ,b ) .  Appendix  1 
gives  an  existence  theorem  for  single-loop  design.  The  first  part  of  the  design 
(see  Appendix  A3)  gives  bounds  on  the  nominal  loop  transmission  which  is 
guu/Puuo  of  w^ere  puuo  the  'nominal'  associated  with  a  nominal  moeft. 
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These  bounds  must  be  satisfied  in  order  that  a  specific  system  transfer 
function  tuy  satisfy  (1).  Here  guu/Puuo  Is  used  for  n  tuy  (v  =  l,...,n) 
functions.  It  is  proven  in  A3,  that  a  9uu/puu0  can  be  found  which  satisfies 
the  conditions  for  all  n  t  functions. 

For  example,  consider  t  ^  at  u  =  u>^  and  suppose  ^ui  (cu-j )  =  .9,  Buj(u». )  = 
1.1  in  (1).  We  could  split  this  range  [.9,  1.1]  into  say  [.95,  1.05]  for  t  ^ 
and  .05  for  T^d^  in  (4),  using  d  -  of  (5)  for  d  -j.  The  technique  in  A3 
or  better  (Horowitz  and  Si di  1972),  is  then  used  to  find  a  bound  on  9uu(jw-|). 
Here,  we  note  a  tough  constraint.  Sooner  or  later  in  u,  |guu(jw)|  must  become 
very  small  with  1  +  9UU/PUU  -*■  1  and  then  in  (4a) 


uv 


q  f  -  d 
auu  uv  uv 


uu 


(10) 


and  in  (7),  ^  ->  the  numerator  of  its  right  side  divided  by  ?uu>  Now  (4a,  5,  6) 
in  general  require  that 

^uvLax  2lTduvduvel  Hi) 

But  ^uv Lax  =  Buv  and  at  high  frequencies 


I Tduvc 


uve 1 


SUPy 

M  ih 


!PuillBivl 

IQ 


To  see  what  this  leads  to  take,  for  example,  n  =  2  so  that  the  above  applied 
to  v  =  1 ,  u  =  1 ,2  gives 


requiring 


1  > 


4 IP1 2P21 

pnp22^ 


as  a)  -*■  ® 


02) 


Thus,  a  constraint  on  P  is 

(P2a):  3  V  ^  for  “  >  %’  1P11P22'  >  4  lP12P2‘l  l^meM* 


03) 


It  is  known  that  as  s  -*■  ®, 


U 


hi. 

,eij 


so  the  above  becomes 


I k1 1  k22  I  >  4Jk12k21 


ell+e22  e12+e21 

U>  OJ 


If  the  uncertainties  in  the  k.^  are  independent  and  e^  +  e22  =  e12  +  e2-| 
this  becomes 


kllmink22min  >  4k12maxk21max' 


04) 


There  is  an  important  problem  class  for  which  the  inequality  is  less 
harsh.  This  is  the  "basically  noninteracting"  class,  where  one  ideally  desires 
t.  •  3  0  for  i  f  j,  but  because  of  uncertainty  accepts  A..  =  0,  |t..|  <  B.. 

1 J  *  J  *  J  '  J 

for  i  f  j,  in  (1).  Also,  one  doesn't  care  if  t.j(i  t  j)  is  nonminimum-phase. 
Condition  (6)  then  applies  only  to  u  =  v.  The  f  (u  f  v)  are  set  equal  to 
zero  and  (13)  becomes 

3  V  ^  '  lPllP22l  >  2  lP12P21  I  meM,  id  >  o)^.  (15) 
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It  is  desirable  to  ease  inequality  (13)  in  the  general  case.  Note  that  (6) 
can  be  satisfied  over  any  finite  w  range  by  making  |1  +  9UU/PUU!  large  enough. 

Thus,  as  previously  indicated,  one  can  split  the  [AUV,BUV]  tolerance  so  that 
l*J  >  lTduvl  lduvel’^m£W’  e-g*  assi9n  lTuvf  e  [E  -  e,  E  +  e]  with 

E  =  ^Auv  +  Buv>/2’  2e  <  BUV  '  Auv  and  the  balance  (Buv  "  Auy  -  2e)/2  is  assigned  to 
Tduvduv  ^a).  But  |1  +  9ul|/£uu!  must  then  be  made  large  enough  to  satisfy 
the  resulting  requirements,  and  it  can  for  any  finite  u>  range.  The  trouble  is 
guu  IIUJSt  a^owe<^  to  zero  as  m  +•«,  leading  to  (13),  etc.,  if  we 
insist  on  (6).  We  could  ignore  (6)  at  large  u>,  say  for  w  >  u)u,  with  Uu  as 
large  as  desired  but  finite,  letting  !tuv |  «  I Tduv I  lduve I  for  “  >  <V  Then 
for  a)  >  (o^,  (11)  is  replaced  by  the  weaker 


It,  Imax  >  It.  d  I 
uv1  1  duv  uve1 

and  for  n  =  2,  (13)  is  then  replaced  by 

(^2b^:  3  “h*  ^  ^or  u  >  uh’  IB1 lB22 I  >  I 2^21 ^ me^ 


(16) 


(17a) 


An  important  question  is  whether  (17a)  is  an  inherent  basic  constraint  in  the 
presence  of  uncertainty,  no  matter  what  design  technique  is  used,  or  is  due 
only  to  this  specific  design  technique.  The  methods  suggested  in  (Rosenbrock 
1974,  Owens  1978)  to  achieve  diagonal  dominance,  may  be  helpful  in  satisfying 
(17a),  but  they  would  have  to  be  extended  to  uncertain  plants.  Note  that  in 
Rosenbrock  1974,  Owens  1978),  diagonal  dominance  is  desired  V  w  e  [0,«>), 
whereas  in  ^jj)  ^  re9uirecl  only  for 

For  the  analog  of  (17a)  at  n  =  3,  it  is  found  that  diagonal  row  dominance 
of  P’1  for  a)  >  ,  is  a  sufficient  condition.  The  necessary  condition  can  be 

written  as 
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(17b) 


•Sr  t 


3  “H*  S  for  “  '  “H  *  lP1JPJl‘  and 

IPnp33l  >  ^  f P1 2P23 ^  +  1P13P22 1^1  P22P31 1  *  lP21P32^ 
which  can  be  written  as, 

^ P1 1P22P33 1  >  I P11P23P32 1  +  iP12P21P33l  +  I P1 2P23P31 1 

+  |P13p22P3il  +  IP1 3P21P32 ^  for  u>  >  u)H  (17c) 

The  latter  has  the  following  intepretation.  Array  the  matrix  P  ^  in  the  usual 
manner,  but  twice -one  under  the  other  as  in  Fig.  3a.  Then  the  terms  on  the 
right  side  of  (17c)  consist  of  the  products  of  the  entries  crossed  by  the 
dashed  lines. 

However,  if  is  so  used,  it  is  no  longer  possible  to  use  Rouche's 
theorem  and  thereby  prove  each  t^  is  minimum-phase.  But  we  can  still  design 
so  that  the  nominal  t^.  are  minimum-phase  and  we  know  from  (6)  that  t.j(jw)  f  0 
for  u>e[0,wH].  Therefore,  from  the  continuity  of  the  zeros  of  t.. .  with  respect 
to  the  parameters  of  the  system,  if  t.,  has  any  right  half-plane  zeros,  they 

'  J 

must  enter  the  right  half-plane  as  shown  in  Fig.  3b.  It  is  unlikely  that  such 
a  zero  which  must  migrate  all  the  way  up  to  ju>^,  should  move  back  into  the 
significant  control  bandwith  region  A.  The  point  is  that  if  right  half-plane 
zeros  are  "far-off",  they  have  little  effect  and  the  system  is  "dominantly" 
minimum-phase. 
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Rouche's  theorem  can  still  be  used  if  we  can  guarantee  that  (6)  is 

satisfied  for  a  semicircle  consisting  of  the  segment  [-ju>H,ju)H]  and  the  right 

half-plane  half- circumference  of  the  circle  of  radius  u>H,  centered  at  the 

origin.  Then,  there  are  definitely  no  right  half-plane  zeros  of  t. .  in  this 

^  J 

half-circle,  and  the  system  is  "dominantly"  minimum-phase  This  is  quite 
practical  in  the  design  technique  of  (Horowitz  and  Si di  1972),  discussed  in  A3. 

3*1  Modification  of  mapping  $ 

Note  that  for  the  "dominantly  minimum-phase"  and  the  "basically  noninter¬ 
acting"  cases,  the  application  of  Schauder's  theorem  in  (2.1),  Eqs.  (7-9),  etc., 
needs  modification,  because  nonmini  mum- phase  tuy(jui)  cannot  be  uniquely 
derived  from  |t  (jw) | .  Redefine  h  e  Huy  of  2.1  to  consist  of  an  ordered 
pair:  h(oj)  as  before  and  q(ui),  the  imaginary  part  of  huy(jaj)  with 
h  =  | huv( joj)  | ;  h  e  Huy  the  same  as  before  but  q(w)  e  C  [0,°°)  with 
0  <  |q(to)  |  <^h(w).  Constraints  2,3  in  2.1  on  h(w)  also  apply  to  q(o)).  Let 
(HQ)uVCC2  [0,®)  denote  the  set  {(h(u>),  q(u>))}  with  ||(h,q)|j  =  ||h||  +  ||q||. 
Obviously,  (HQ ) uV  is  compact  and  convex  in  C2  [0,®).  The  extension  to  the 
n2  product  set  is  straightforward. 

The  mappings  4>uv  in  (7)  are  redefined.  Each  ^uy  is  a  pair  of  mappings, 
one  the  absolute  value  as  before,  the  second  the  imaginary  part  with  the 
absolute  bars  on  the  right  removed.  On  the  right  side  of  (7),  B(h.y(u>))  is 
replaced  by  r.y(u>)  +  jq . y(o>) ,  with  h?y  =  r.y  +  qiy,  (h.y,  qiv)  e  (HQ)iy. 

It  is  necessary  to  prove  that  i>  maps  each  element  of  (HQ)UV  into  itself. 
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The  proof  follows  immediately  from  that  for  the  minimum-phase  case  --  this 
is  obvious  from  (6),  the  definition  of  duve  in  (5),  and  Appendices  1,2.  The 
proof  that  <P  is  continuous  is  straightforward.  Accordingly,  the  Schauder 
conditions  are  satisfied  and  there  exists  a  fixed  point  which  satisfies  the 
specifications.  Such  specifications,  by  themselves,  would  not  be  good  ones 
because  they  permit  highly  nonminimum-phase  tuv(s).  However,  they  are 
satisfactory  if  it  is  known  from  other  sources  that  t  is  "dominantly 
mini mum- phase". 

3.2  Additional  Constraints  on  P 

Constraints  Al(l)-(3)  in  the  Appendix,  must  be  applied  to  the  1|P  , 

since  in  Fig.  (2)  p  =  =  p  of  Appendix.  Al.l  requires  that  there  be 

u  uu  1  A 

no  change  in  the  excess  of  poles  over  zeros  of  p —  =  where  A  =  det.  P 

uu  fluu 

and  Auu  its  uuth  minor,  as  m  ranges  over  Also,  that  for  at  least  one 
mejl,  denoted  by  muQ,  Pu(j  has  all  its  poles  and  zeros  in  the  interior  of  the 
left  half-plane.  The  mUQ  can  be  different  for  each  u. 

A1.2  requires  that  1/PUU  is  minimum-phase V  me/t,  and  its  zeros  do  not 
get  arbitrarily  close  to  the  jw  axis.  Since  1/PUU  =  A/Au(J,  this  means  A  must 
have  no  right  half-plane  zeros.  Hence  the  P..  in  general  have  no  right 

*  J 

half-plane  poles.  (For  those  who  wish  it,  P  is  restricted  to  be  controllable 
and  observable  V  mefl,  but  these  concepts  are  unnecessary  if  P  is  properly 
formulated  in  terms  of  physical  uncertain  parameters  (Horowitz  and  Shaked 
1975)).  Since  the  p..  in  P  =  [p..]  are  finite  rational  functions,  the  latter 

I  vJ  I  J 

part  of  A1.2  is  automatically  satisfied. 


-94- 


*>.'  * 


7* 


A1.3  for  n  =  2  is  the  same  as  (17),  which  shows  that  (17)  is  a 
fundamental  condition  for  linear  time-invariant  design,  not  an  "extra" 
condition  due  to  our  design  technique,  at  least  for  n  =  2.  However,  (13) 
is  an  "extra"  condition.  Note,  the  extension  of  single-loop  design  to 
disappearing  poles  and  zeros  in  A6  may  perhaps  permit  disappearing  poles 
and  zeros  in  the  mio  plant  functions. 


4.  OTHER  DESIGN  EQUATIONS 


The  previous  design  equations  constitute  only  one  of  many 
design  techniques  derivable  from  Schauder's  fixed  point  theorem.  Only 
two  more  will  be  briefly  mentioned  here. 

Both  are  based  on  the  use  of  a  nominal  diagonal  loop  trans¬ 
mission  matrix.  The  design  obligations  on  the  loop  transmission 


elements  are  then  independent  of  the  way  the  plant  input  and  output 
terminals  are  numbered.  If  G  is  made  diagonal,  such  numbering  is 
important  and  after  one  arbitrarily  numbers  the  plant  input  terminals, 
he  should  try  to  number  the  outputs  such  that  the  main  effect  of  in¬ 
put  i  is  on  output  i.  Manipulation  of  (2)  somewhat  differently  from 
Sec.  2,  gives 


fn*n/6n  *  m  11  1]/  11 

ii  “  i  + 

t,.  -  V22'«22  +  W  V2itll/*22>  6tC' 

1+V622 


Lo  '  PoG 

Po  =  t  V 


(18) 
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where  V  =  [v^]  =  I-.^(P)”^,  Pq  is  the  'nominal' plant  matrix  and  there¬ 
fore  fixed,  P  is  the  general  uncertain  plant  matrix,  6.^  =  1  -  v^. 

The  «...  are  the  nominal  elements  of  the  loop  transmission  matrix  L. 
n  r 

Eqs.  (18)  lend  themselves  to  single-loop  design  and  use  of  Schauder's 
theorem,  precisely  as  did  Eqs.  (4). 

Another  interesting  set  of  design  equations  is  obtained  by  de¬ 
signing  to  control  the  changes  in  t^,  rather  than  t.^  directly.  Let 

T  =  [t..  ]  be  the  'nominal'  system  transfer  matrix  and  T  =  [t..]  the 
o  1JO  1J 

actual  which  is  uncertain,  AT  =(At..]  =  T  -  T  .  Then  it  can  be  shown 

U  o 

that 


AT  =  (I+L)'1  VT,  V  =  1-2 0P_1 


(19) 


where  2  ,  2  are  likewise  the  'nominal'  and  uncertain  plant  transfer 
o 

matrices,  and  L  =  I  G  =  [A..]  is  the  nominal  loop  transmission  matrix, 
o  ij  - 

If  L  is  taken  diagonally,  the  result  is  (n  +  2  for  simplicity) 


At 


1.11  T  12  21 


11 


1  +  l 


At 


1112  12  22 


11 


12 


1  +  i 


(20) 


11 


and  similar  obvious  ones  for  At^,  At^. 

The  design  problem  is  now  completely  one  of  disturbance 
attenuation,  with  the  disturbances  =  v^t^  +  v^t^,  etc.,  whose 
range  is  known.  Schauder's  theorem  is  applicable  in  the  same  manner 
as  before.  Note  that  V  represents  the  'normalized'  plant  variation 
matrix.  Eqs.  (20)  appear  to  be  much  simpler  to  use  for  design  (once 
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the  At.,  tolerances  are  formulated)  than  (4),  and  their  use  needs  to 

*  J 

be  intensively  researched.  However,  both  for  (18)  and  (20)  the  con¬ 
straints  considered  in  3.,  leading  to  (11-15)  must  be  found,  and 
these  may  possibly  be  tougher  than  before.  Also,  both  a  nominal  £ 
and  T  must  be  chosen^ which  is  not  goodj because  the  optimum  pairing 
is  not  apriori  known.  However,  the  analogs  of  (14,17)  may  be  more  lenient. 

4.1  Bandwidth  Minimization 

An  important  criterion  for  comparison  of  design  techniques  is 
their  "cost  of  feedback,"  which  we  take  as  the  bandwidths  of  the 
loop  transmission  functions— because  they  determine  the  system 
sensitivity  to  sensor  noise.  Obviously,  quantitative  synthesis 

techniques  must  first  be  invented  before  one  can  turn  to  their  op¬ 
timization  (for  without  such  quantitative  techniques  comparison  is 
possible  at  best,  by  analysis  after  a  specific  numerical  design  has 
been  made).  This  approach  via  Schauder's  theorem  promises  to  generate 
a  variety  of  such  techniques,  and  the  next  step  will  be  optimization. 


5.  DESIGN  EXAMPLE 


The  2x2  plant  elements  are  p..  =  k../0+sA..)  with  correlated 

1J  i J  * J 

uncertainties,  giving  a  total  of  9  parameter  sets  in  Table  1.  The 
design  was  performed  to  handle  the  convex  combination  generated  by 
these  9  sets  (Figure  6). 
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TABLE  1 


* 


No. 

kn 

^22 

k12 

k21 

5li 

^22 

A12 

^21 

1. 

i 

2 

.5 

1 

1. 

2 

2 

3 

2. 

i 

2 

.5 

1 

.5 

1 

1 

2 

3. 

i 

2 

.5 

1 

.2 

.4 

.5 

1 

4. 

4 

5 

1 

2 

1. 

2 

2 

3 

5. 

4 

5 

1 

2 

.5 

1 

1 

2 

6. 

4 

5 

1 

2 

.2 

.4 

.5 

1 

7. 

10 

8 

2 

4 

1. 

2 

2 

2 

8. 

10 

8 

2 

4 

.5 

1 

1 

2 

9. 

10 

8 

2 

4 

.2 

.4 

.5 

1 

A  " 

basically 

noninteracting" 

system  is 

desired. 

with 

the  off-diagonal 

transmissions  specified  in  the  (^domain  |t^(jw)|,  |t2^{jw)|<0.1  Vw.  The 
diagonal  t^,  bounds  are  identical  and  were  originally  in  the  time- 
domain  in  the  form  of  tolerances  on  the  unit  step  response  shown  in 
Fig.  4a,  b  (which  also  shows  the  design  results  for  those  of  the  9  cases 
which  were  reasonably  distinguishable).  These  time-domain  bounds 
were  translated  into  the  "equivalent"  bounds  on  |t.,.(ju))|  shown  in 
Fig.  5  (Horowitz  and  Sidi  1972,  Krishnan  and  Cruickshank  1977). 

Familiarity  with  quantitative  single-loop  design  is  assumed 
here  .  One  can  do  a  problem  of  this  complexity  by  hand.  The 

sets  {p..  (jw)},  called  the  plant  templates,  are  obtained  on  the 

1ie  -i 

-1  A  xF*  A 

Nichols  chart.  Some  of  these  templates  of  =  - — , =  « — 

11  '22  ’ll 
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are  shown  In  Fig.  6  at  various  u>  values.  The  larger  the  template. 


the  greater  uncertainty  at  that  ui  value.  The  tolerances  on  t  of 

(4a)  and  Fig.  5  were  divided  between  t  and  t  d  as  discussed  In 
'  3  uu  duu  uu 

Sec.  2.  Each  of  these,  in  conjunction  with  the  templates,  leads  to 

9.... 


bounds  on  the  nominal  loop  transmission  i 


uuo  2 


Some  of  these 


bounds  on  ,  due  to  t  ,  are  shown  as  solid  lines  in  Fig.  7,  i.e. , 
no  11 

it  is  necessary  for  to  lie  above  the  indicated  boundary.  The 
tolerances  on  t.  d  lead  to  the  dashed  line  bounds  on  JL,  .  No 
attempt  was  made  to  optimize  the  division  of  the  tolerances  between 
T-j .j  and  Tdl-|dyi-  The  composite  bound  on  £^o  must  satisfy  both. 

The  &n0(jw)  chosen  is  also  shown  in  Fig.  7.  There  was  no  attempt 
made  to  optimize  the  JL.q;  the  design  was  made  by  hand  quickly,  so 
the  are  larger  than  need  be,  with  the  tolerances  therefore 

satisfied  better  than  necessary— as  seen  in  Figs.  4a,  b.  Optimal 
£^0(jw)  would  lie  on  their  boundariesat  each  o> ,  so  in  this 
example  there  is  considerable  overdesign. 

Here  we  took 


o  10  (H-.007S) _ 

llo  "  p_.  911  ~  s  (1+.025s)M+s  , 

l  400  (400)2 
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wi  th 


with 


Ao  .75  (1+3. 66s) 
f22o  "  (l+s)(l+3s)  *’ 


22o 


llo 


’22 


=  9  (1+.  02s) 


O+.ls) 


1+  s 
.  150 


(150) 


2\ 


Ao  1.5  (1+3. 66s)  . 

fllo  '  (l+3s)(l+2s) 

The  requirements  on  f^,  f ^  <fi2  =  f21  =  9 12  =  g21  =  0  here^  were 
found  using  single-loop  design  technique  [15]  as  briefly  explained 
here  in  A4,  and 

f  =  — 3 —  f  =  — 1 _ 

11  1  +  .5s  *  r22  1  +  .33s 

were  found  satisfactory.  The  system  was  simulated  on  the  digital  com¬ 
puter  with  the  results  shown  in  Figs.  4a,  b.  The  t^»  t^^  tolerances 
were  easily  satisfied  by  the  design. 

While  this  is  not  a  very  challenging  example  of  the  design  techni¬ 
que,  nevertheless  the  uncertainty  is  very  large  and  one  should  consider 
how  quick,  simple  and  straightforward  was  the  design  procedure,  and 
also  consider  what  alternatives  are  offered  in  the  mio  literature. 

There  are  no  other  techniques  available  for  systematic  design  to 
specifications  in  the  presence  of  significant  uncertainty  which 
guarantee  design  convergence  and  attainment  of  design  tolerances. 


-100- 


p 


Whatever  present  popular  technique  is  used,  it  would  be  necessary 
to  cut  and  try  and  endeavor  to  understand  the  relations  between  the 
cutting  and  the  results  as  one  continued  to  cut  and  try,  because 
these  techniques  have  no  provision  for  significant  uncertainty.  In 
the  above  design,  one  sweep  was  known  to  be  sufficient  because  the 
plant  and  the  design  tolerances  (w-domain)  satisfied  constraints, 

PI  etc. 


5.  EXTENSION  TO  NONLINEAR  UNCERTAIN  MIO  PLANTS 

Once  there  is  a  quantitative  design  technique  for  linear  time 
invariant  mio  uncertain  plants,  it  appears  at  least  conceptually 
possible  to  extend  it  to  a  significant  class  of  nonlinear,  even 
nonlinear  time-varying,  uncertain  mio  plants.  The  procedure  is  a 
generalization  of  that  used  (based  also  on  Schauder's  theorem) 
in  (Horowitz  1976)  for  single  loop  uncertain  nonlinear  systems. 

The  key  feature  is  the  replacement  of  the  nonlinear  plant  matrix 
set  (a  set  because  of  the  uncertainty),  by  a  linear  time  invariant 
plant  set  which  is  precisely  equivalent  to  the  original  nonlinear 
set,  with  respect  to  the  acceptable  system  output  set.  The  pro¬ 
cedure  is  briefly  presented  for  the  case  where  one  wants  the  system 

with  nonlinear  uncertain  plant  to  behave  like  a  linear  time-invariant 
system  for  a  specified  class  of  command  input  sets. 
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It  is  essential  that  the  command  input  sets  represent  a  good 
sampling  of  how  the  system  will  actually  be  used.  For  example,  suppose 
n  =  3  and  in  actual  use  r^,  r^  always  exist  simultaneously  (with 
r^  a  0),  and  r3  appears  by  itself  (with  r^  =  r„  =  0).  Say  there  are 
ten  typical  r^(t)  inputs  and  for  each  typical  r^(t)  there  are  five 
typical  r2(t).  This  makes  a  subtotal  of  50  input  sets,  to  which  is 
added  the  number  of  typical  r3(t)  say  10,  giving  a  class  R  =  { r)  of 
60  sets,  of  which  50  have  the  form  r  =  (r^,  r 2>  0)  and  r  =  (0,  0,  r3) 
for  the  balance.  Choose  r^eR.  The  family  of  acceptable  outputs  for 
this  input,  is  known  from  the  tolerances  on  t..,  giving  for  that  one 
input  vector  a  family -H  -  {h},  h  =  (h^ ,  h2»  h3).  The  mio  plant  is  re¬ 
presented  by  a  family  (because  of  parameter  uncertainty)  W  of  nonlinear 
differential  mappings 

U  -  M,  W  =  (w] ,w2,w3)^c1  =  w1(x2,x2,x3>m),  .  .  .  ,  c3  = 
w3(xl ,x2,x3,m),  where  the  x*  are  the  plant  inputs,  c^  the  plant  outputs, 
and  m  is  the  plant  parameter  vector  meM. 

Take  a  sample  acceptable  outputtriple  h  =  (h^,  h2>  h3)  and  find 
the  corresponding  plant  inputs  at  some  specific  meM  (or  in  other  words, 
pick  a  Weto)  and  let  c^  =  h^  and  solve  the  nonlinear  equations  backwards, 
giving  the  input  set  (x^ ,  x2>  x3).  Take  the  Laplace  transforms  x^(s) 
of  X.,  fij(s)  of  hj  giving  the  vectors  x[s]  =  (x^s),  x2(s),  x3 ( s ) ) 

fits]  =  (fi^s),  ....  fi3(s)).  Repeat  for  other  h  samples  in  the 
acceptable  output  set  H,  giving  two  paired  families  of  fc[s],  h . [s] . 


«».-■> '  w 


Select  any  combination  of  three  x[$],  forming  a  3  x  3  matrix  X  and 

A 

corresponding  paired  combination  of  three  n[s],  forming  the  matrix  H. 

Set  H  =  PX  and  solve  for  P  =  H(X)"1.  P  is  the  linear-time-invariant 
equivalent  of  the  specific  WeW  picked,  with  respect  to  the  specific 
trio  of  acceptable  output  vectors  picked.  Repeat  over  different 
trios.  Repeat  the  entire  operation  over  different  wefc',  giving  a 
class  Pj  *  {P}  which  is  the  linear-time-invariant  equivalent  of  the 
W  family,  with  respect  to  the  class  of  acceptable  outputs  H  for  in¬ 
put  vector  ?i .  Repeat  the  entire  operation  for  ry  .  .  .  ,  rg0  giving 
{?.}  =  P^  ^  ,  which  is  the  linear  time  equivalent  for  the  nonlinear 
W,  with  respect  to  the  tribe  of  60  families  of  acceptable  output  sets. 
The  equivalence  is  exact  if  the  conditions  for  application  of 
Schauder's  theorem  are  satisfied.  We  now  have  a  linear  time-invariant 
uncertain  mio  problem,  which  let  us  presume  we  can  solve.  If  and  only 
if  we  can  guarantee  the  solution  of  the  latter,  then  the  same  compen¬ 
sation  functions  will  work  for  the  original  nonlinear  uncertain  mio 
plant.  Hence  the  importance  of  quantitative  linear  time  invariant 
design  techniques  (over  and  above  their  intrinsic  importance)— for  they 
enable  the  precise  solution  of  nonlinear  uncertainty  problems. 

The  design  effort  in  the  above  appears  to  be  enormous  but  it  is 
conceptually  straightforward  and  easy.  An  ordinary  control  engineer 
can  implement  it  and  the  digital  computer  is,  of  course,  an  essential 
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tool.  Conceptually  too,  it  appears  possible  to  extend  the  method  to 
obtain  nonlinear  relations  between  inputs  and  outputs  within  specified 
bounds,  despite  large  plant  uncertainty,  even  nonlinear  time-varying, 
as  can  be  done  for  the  single  input-output  case.  The  prospect 

is  fascinating.  Imagine  being  able  to  work  with  the  actual  nonlinear 
equations  of  a  jet  engine,  or  a  chemical  process,  etc.,  include  un¬ 
certainties  in  the  modelling,  even  uncertainty  in  system  order  (see 
Appendix),  and  designing  to  achieve  outputs  within  specified  tolerances 
over  the  given  range  of  uncertainty. 


* 


L 


K  n  )» 


6.  DISTURBANCE  ATTENUATION 

Let  x  in  Fig.  1  be  a  n x 1  matrix  of  disturbances.  The  resulting 
system  output  (with  r  =  0)  is  c  =  (I  +  PG)  ^  Px  -  Zx,  Z  =  [z.jg]»  nxn 

disturbance  response  matrix.  Bounds  on  Z  are  given  in  the  form 


z  (jw)|  <  buv(u>)  .  V  meM 


(21) 


Rewrite  c  =  Zx  in  the  form  (P"1  +  G)c  =  x.  Let  x.  i  0  only  for  i  =  v,  so 


ci  =  zivV  and 


M 

r  ,  *  .u  _  /0,  u  f  V\ 

£  ^ui  ^ui  ziv  -  v  1,  u  =  v  ’ 


i=l 


<puu*9UU)2Uv  '  •  Z  tpui  + 


1/U 


5v  ^  ^Pui  +  9ui ^zi\ 

1£h _ 


'uv 


p  (i  +  ) 

uuv 1  +  p  ) 
uu 


(22) 
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xuve(u))  '  sup  l  I"511? — — lbiv(a>)  (23) 

m  i/u 

The  9u1-(w)  (i^u)  can  be  chosen  to  minimize  xuve(i»)>  but  for  simplicity  we  shall 
assume  them  zero.  From  (22,23) 


uv 


(u))|  < 


,VPuu  +  xuve| 

(i 

uu 


(24) 


If  1/P  satisfies  the  constraints  listed,  then  it  is  obviously  possible 

to  guarantee  |zuv(u>)|  <  any  finite  number,  no  matter  how  small,  at  any  finite 

(jo.  Also  it  can  be  made  zero  at  a  finite  number  of  u>  values  by  assigning  poles 

to  guu  at  these  values.  Assume  that  guu  can  be  chosen  to  satisfy  (21) 

V  a)  e  [0,”).  Then  one  can  set  up  the  conditions  for  Schauder’s  theorem, 

precisely  as  was  done  in  2.1.  The  set  buv(u>)  must  have  been  formulated  such 

that  B(n  ),  the  n  product  set  of  the  buy(aj),  is  compact  convex  in  C(n  )v 
o 

analogous  to  H(n  )  in  2.1.  The  analog  of  4>  in  (7)  must  be  formulated 

with  the  modification  of  Sec.  3.1,  inasmuch  as  we  do  not  care  if  the  zuv(s)  are 

nonmini mum- phase. 

Conditions  analogous  to  (12-17)  for  n  =  2,  are  obtained  as  follows. 

As  w  -*■  °°,  9UU/PUU  -*■  0  so  in  (24),  the  right  side  -►  its  numerator.  But 
|zuv(jw)|  <  buv(u)  of  (21).  Let  u  =  1,  v  =  2  and  then  u  =  2,  v  *  1  and 
obtain  the  necessary  condition  (for  g12  s  92]  = 


As  a)  -*  °°  , 


p12p21  <  pnp22  ’  ^ 


(25) 


f  » 

I 


i 


similar  to  (17)  but  here  only  at  °°,  because  there  is  no  concern  with  the 
minimum-phase  property.  Setting  u  =  v  =  1,  and  then  u  =  v  =  2  in  (24),  we 
get  the  conditions 


As  a) 


'll 


11 


ip,i  - 


p12p21 

p22 


b22  >  ^22- 


p12p21 

p11 


(26) 


But  in  reality  as  u  -*  °°,  c  -*■  Px  so  Z  -*■  P  and  -*•  P-|i»  z22  P22*  ^ence, 

assignment  of  b^  (as  w  ■>  °°)  to  satisfy  (25)  is  no  obstacle,  because  the  buv(u>) 
are  upper  bounds  on  the  |zuv(joj)|. 
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